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ABSTRACT 

An  analytical  investigation  is  undertaken  into 
the  application  of  the  upper  bound  theorem  of  limit 
design  for  use  in  the  design  of  ship  transverse  web 
frame  structures.   The  revised  simplex  method  is  used 
to  solve  the  upper  bound  theorem  relationships  for 
required  values  of  plastic  limit  moment  for  structurally 
safe  designs.   Bracket  and  beam  overlap  influence  is 
idealized  by  assuming  rigid  extensions  in  the  overlap 
and  bracket  regions.   Comparisons  are  made  with  existing 
web  frames  designed  by  elastic  methods. 
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NOTATION 


A  Total  beam  cross  sectional  area. 

A„  Beam  Flange  Area. 
F 

Am  Cross  sectional  area  of  beam  web 
T 

and  top  flanges. 

A  Coefficient  matrix  of  constraint 

equation. 

a  Distance  between  centroids. 

aij  Constants  of  constraint  function. 

B  Vector  of  the  basis. 

b  Right  hand  side  of  constraint  relation. 

b-  Constant  of  constraint  function. 

c^f  c  •  ,  c^          Length  of  beam  overlap. 

c  Coefficient  matrix  of  objective  function, 

d.  Bracket  leg  length. 

f  Objective  function. 

fT,  f               Lengths  of  rigid  sections  due  to  beam 

overlap. 

g  Beam  element  weight  per  unit  length. 

h  Beam  depth. 

L.,  1.              Beam  element  length. 


M 


a,  Mfc,  M. ,  M0      Plastic  limit  moment. 


P.  Vector  associated  with  extreme  point. 

-3 

Sy,  s  Overlapping  web  depths. 

S2  Cross-tie  height  above  base  line. 

S5  Side  tank  width. 
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NOTATION  (cont.) 

t£  Rigid  extension  due  to  bracket 

influence. 

u  Normal  velocity. 

w  In  line  velocity. 

x.  Variable  of  linear  programming 

J  problem. 

z.:  Iteration  solution  value  of  objective 

function. 

z  Plastic  section  modulus. 
P 

tt_  Multiplier  vector. 

p  Material  density. 

av  Material  yield  stress. 
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1.   INTRODUCTION 

1.1   Preface 

The  design  problem  of  transverse  web  frame 
structures  in  ships  is  receiving  an  increasing  amount  of 
interest  with  the  advent  of  mammoth  tankers  and  bulk 
carriers.   Transverse  framing  accounts  for  a  sizable 
portion  of  the  total  structural  weight  of  such  ships.   A 
substantial  weight  savings  can  result  from  the  proper 
design  and  optimization  of  the  repeating  transverse 
frames  in  the  prismatic  mid  section  of  these  large  vessels. 
Great  effort  is  usually  made  in  the  design  phase  to  make 
these  structures  as  efficient  as  possible.   Modern  day 
design  efforts  have  been  concentrated  in  two  methods  of 
elastic  design. 

The  finite  element  method  is  the  most  popular  of 
the  modern  elastic  ship  design  methods.   A  number  of 
review  papers  [l  -  4]  have  discussed  the  application  of 
this  technique  to  the  design  of  ship  structures  in  general. 
Smith  &  Woodhead  [l]  discussed  specific  application  of  the 
method  to  web  frames. 

The  second  elastic  design  method  uses  beam  theory 
and  was  first  discussed  in  depth  by  Lund  [5].   Others  have 
also  written  on  the  topic  [4  -  14]. 

Both  of  the  elastic  design  methods  involve 
systems  of  nonlinear  equations.   Design  optimization  can 
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become  complex  and  expensive  in  terms  of  computation  time. 
The  nonlinearities  require  the  use  of  optimization 
techniques  such  as  nonlinear  programming,  Powell's  direct 
search  method  or  the  gradient  method. 

The  elastic  design  approach  and  the  nonlinear 
optimization  methods  combine  to  form  a  laborous  design 
task  with  a  trusted  but  possibly  expensive  and  inefficient 
product. 

The  purpose  of  this  study  is  to  examine  the  use 
of  another  design  method,  namely  plastic  limit  design, 
in  the  design  of  transverse  web  frames  in  ships. 
1. 2   Historical  Review 

A  search  of  existing  literature  reveals  a  wealth 
of  information  on  the  use  of  plasticity  design  methods  and 
specifically  the  use  of  the  plastic  limit  design  theorems 
dating  back  several  decades  [15  -  22],   The  limit  design 
methods  are  common  place  and  vital  in  nearly  all  sectors 
of  structural  design  except  naval  architecture  and  ocean 
engineering.   In  papers  on  marine  structures  where  the  use 
of  limit  design  is  suggested  [23  -  25],  it  is  often  from 
the  civil  engineering  point  of  view  with  examples  taken 
from  standard  civil  engineering  texts. 

Jones  [24,  25]  uses  the  limit  theorems  in  the 
analysis  of  tanker  web  frames  in  two  separate  papers.   To 
date,  his  are  the  only  papers  known  to  discuss  in  detail 
the  use  of  the  limit  theorems  on  marine  structures  loaded 
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by  various  conditions  of  hydrostatic  loading. 
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2.   METHOD  OF  CALCULATION 

2. 1   Limit  Design  Upper  Bound  Theorem 

In  structures  composed  of  an  elastic,  perfectly 
plastic  material,  one  of  the  primary  functions  of  the 
structure  is  that  it  should  carry  a  given  set  of  loads 
while  in  equilibrium  without  flow  occurring  in  the 
structure.   The  structural  design  problem  is  referred  to 
as  the  limit  design  problem  if  this  is  the  only  function 
of  the  structure  [22].   An  alternative  formulation  of  the 
limit  design  problem  states  that  if  under  all  possible 
failure  modes  the  internal  energy  dissipated  during 
kinematically  admissible,  small  virtual  displacements 
exceeds  the  external  work  rate  associated  with  these 
displacements  then  the  structure  is  safe  [7].   This 
latter  statement  can  be  used  to  formulate  an  upper  bound 
on  the  system  plastic  moments  as  a  programming  problem. 

The  limit  design  upper  bound  theorem  is  the 
primary  tool  to  this  analysis.   The  desired  output  of  the 
limit  design  problem  is  the  required  weight  (i.e.  cost)  of 
individual  structural  elements  and  the  total  system. 

The  output  of  the  limit  design  process  is  the 
plastic  moment,  therefore,  it  is  necessary  to  establish 
some  relationship  between  weight  per  unit  length  g  and 
plastic  moment  M0.   Home  [21]  states  that  for  materials 
of  given  yield  stress  and  geometrically  similar  section, 
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g  «  m2'3.      For  Universal  Beam  Sections  the  relationship 
has  been  shown  to  closely  approximate  g  «  M0*6.   More 
generally  g  a   Mn.   Home  also  states  that  the  relative 
plastic  moments  for  absolute  minimum  weight  in  a  given 
structure  are  only  very  slightly  affected  by  the  exact 
value  chosen  for  the  index  "n".   It  will  be  shown  that 
solutions  are  greatly  simplified  by  maintaining  linear 
relationships  between  variables.   It,  therefore,  is 
assumed  that  the  weight  per  unit  length  of  an  element 
is  given  by 

g  =  AL  +  BLM. 
where  A  and  B  are  constants  and  L  is  the  element  length. 
The  total  system  weight  of  the  structure  is  then  given 
by 

G  =  I    A±  L±    +  I    Bi  L.  Mi 
The  methods  of  plastic  limit  design  cannot 
utilize  the  principle  of  superposition  of  stresses  and 
deflections  due  to  individual  loads,  as  is  possible  in 
design  where  the  structure  remains  in  the  elastic  range 
[15 J.   A  paradox  of  the  limit  theorem,  however,  is  that 
when  a  collapse  mode  of  a  structure  involves  only  part 
of  the  structure,  the  collapse  load  of  the  plastic  limit 
moment  can  be  expressed  independently  of  the  characteristics 
of  the  remainder  of  the  structure  [18].   In  other  words, 
in  formulating  the  limit  design  problem,  only  those  loads 
associated  with  a  given  collapse  mechanism  need  be 
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considered  with  that  mechanism.   This  premise  allows  the 
energy  equations  associated  with  each  mechanism  to  be 
constructed  independently. 

Further  assumptions  and  idealizations  are 
required  to  complete  the  formulation  of  the  limit 
design  problem.   Material  behavior  is  idealized  to  be 
elastic,  perfectly  plastic  (i.e.  free  of  strain  hardening) 
Strains  are  assumed  to  be  small.   Jones  has  discussed 
the  importance  of  these  assumptions  [26], 

Finally,  it  is  assumed  that  all  frames  are 
plane  frameworks  of  beams,  in  which  transverse  shear  and 
axial  forces  do  not  influence  the  plastic  yielding  of  the 
beam.   Symonds  and  Neal  [15]  stated  that  shear  forces 
can  be  neglected  except  in  cases  where  the  ratio  of  span 
length  to  beam  depth  is  less  than  approximately  four  to 
one.   Also,  the  effect  of  an  axial  force  on  the  fully 
plastic  moment  can  be  neglected  when  the  force  is  less 
than  about  one  fourth  of  the  compressive  yield  force 
of  the  member. 

To  illustrate  the  general  proceedures  of 
formulating  a  limit  design  problem,  the  simple  civil 
engineering  portal  frame  of  Figure  1,  will  be  examined. 
Application  to  the  more  sophisticated  ship  frames  with 
hydrostatic  loads  is  discusses  in  Reference  [25]  in  the 
appendices  and  in  sections  2.2  to  2.4.   The  frame  is 
restrained  vertically  as  well  as  horizontally  at  its 
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base  and  the  supports  are  clamped  against  rotation. 
The  frame  is  externally  determinate  and  internally 
indeterminate  with  r  =  3  redundancies.   A  discussion  of 
redundancies  may  be  found  in  [27], 

The  frame  has  n  =  8  critical  sections,  or 
locations  where  plastic  hinges  can  form.   Plastic  hinges 
may  form  at  the  ends  of  a  span,  directly  under  a 
concentrated  load  or  within  the  span  of  a  member 
subjected  to  a  distributed  load. 

Only  r  of  the  n  critical  bending  moments  can 
be  defined  independently  of  the  equilibrium  equations 
[24 J.   Therefore,  there  exist  n  -  r  or  5  independent 
relations  between  the  bending  moments.   These  relationships 
are  the  system  equilibrium  equations  and  can  be  obtained 
by  utilizing  the  principle  of  virtual  work.   Each 
equilibrium  equation  is  associated  with  a  kinematically 
admissible  fundamental  mechanism.   Fundamental  or 
elementary  mechanisms  are  classified  by  Symonds  and  Neal 
as  beam,  panel,  sway  and  joint  rotation.   Joint 
rotations  are  also  classified  fictitious.   The  5 
mechanisms  of  the  sample  problem  are  illustrated  in 
Figure  2.   For  all  5  failure  mechanisms  to  safely 
carry  the  prescribed  loads  we  must  have 

M4  -  2M5  +  M6  >  2  (2.1) 

M-l  -  4M2  +  3M3  >  6  (2.2) 

M,  -  M3  +  M?  -  M8  >  6  (2.3) 
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M3  "  M4  -  °  (2.4) 

M6  -  M?  >  0  (2.5) 

where  values  of  6,  angular  velocity  have  been  divided 
from  both  sides  of  the  inequalities.   Minus  signs  denote 
hinge  rotations  in  which  the  internal  fibers  of  the 
hinges  are  in  tension.   This  convention  is  merely  a 
method  of  keeping  track  of  the  direction  of  rotation 
and  is  not  intended  to  imply  negative  internal  energy 
dissipation. 

The  fundamental  sway  mechanism  was  eliminated 
from  consideration  in  this  study.   This  was  not  because 
of  any  assumptions  that  sway  would  not  cause  collapse 
in  ship  web  frames,  but  rather  that  sway  could  not  be 
induced  by  any  of  the  symmetrical  loading  models  used 
by  the  elastic  designs  which  are  used  as  comparison 
baselines  for  this  study. 

The  fundamental  mechanisms  are  not  the  only 
mechanisms  which  must  be  examined.   It  is  possible  for 
linear  combinations  of  the  fundamental  mechanisms  to 
control  the  design  [28],   This  is  true  where  the  combination 
of  mechanisms  allows  a  reduction  or  elimination  of 
certain  plastic  hinges  reducing  the  total  internal 
energy  dissipation  while  the  combined  external  work 
rate  remains  unchanged.   It  is  here  that  the  signs 
of  the  various  terms  in  the  equilibrium  equations  become 
helpful  bookkeeping  aids.   In  combining  the  fundamental 
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mechanisms  of  the  example  the  fictitious  joint  rotations 

are  used  to  derive 

M3  -  2M   +  M?  >  2  (2.6) 

M1    -  4M2  +  3M4  >  6  (2.7) 

Ml  "  M4  +  M7  "  M8  -  6  (2.8) 

also  equations  (2.2),  (2.3)  and  (2.5)  can  be  combined  to 
give 

4M1  -  4M2  +  3M7  -  3Mg  >  24     (2.9) 
or         4MJL  -  4M2  +  3M6  -  3Mg  >  24     (2.10) 
Assume,  now,  that  the  structure  will  be 
developed  with  two  separate  beam  sizes.   The  two  vertical 
elements  will  be  identical  in  cross  section  with  plastic 
moment  M  and  the  horizontal  element  is  allowed  to  differ 

cl 

from  the  verticals  with  plastic  moment  M,  .   Now  equations 
(2.1)  through  (2.3)  and  (2.6)  through  (2.10)  can  be 
written  in  terms  of  two  variables  with  the  sign 
designators  removed  as  follows 

4  Mb  >  2  (2.1a) 

8  M   >  6  (2.2a) 

4  Ma  >  6  (2.3a) 
2  Ma  +  2  Mb  >  2  (2.6a) 

5  Ma  +  3  Mb  >  6  (2.7a) 

2  M   +  2  M,  £  6  (2.8a) 

a      b 

14  M  24  (2.9a) 

a.  ~ 

11  Ma  +  3  Mb  >  24  (2.10a) 

Note  that  the  fictitious  joint  rotations  can  not  physically 
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occur  alone  and  therefore  appear  in  this  final  formulation 

only  in  the  formation  of  the  combined  mechanisms.   The 

minimum  weight  design  problem  is  then  to  minimize 

Min  f  =  2  L   M   +  L,  M, 
a   a    b  b 

or   Min  f  =  8  Ma  +  4  Mfa  (2.11) 

subject  to  expressions  (2.1a)  to  (2.3a)  and  (2.6a)  to 
(2.10a),  where  the  function  (2.11)  is  the  design  minimum 
weight  function. 

It  is  here  that  the  assumption  of  a  linear 
weight  and  moment  relationship  discussed  on  page  4 
becomes  important.   In  the  present  linear  form  the 
problem  of  the  minimum  weight  design  becomes  a  relatively 
simple  linear  programming  problem  in  two  variables,  M 
and  M,  ,  where  the  object  function  is  (2.11)  and  expressions 
(2.1a)  to  (2.3a)  and  (2.6a)  to  (2.10a)  are  the  constraining 
functions. 

This  linear  programming  problem  will  be  solved 

for  M„  and  M,  in  section  3.1. 
a      b 

2. 2   Rigid  Bracket  Region  Determination 

In  the  limit  design  problem  of  very  large  and 
deep  ships,  with  deep  web  frame  structures,  the  effect 
of  web  overlap  and  bracket  influence  becomes  significant. 
It  is  customary  [5,  14]  in  elastic  design  of  ship  frame 
structures  using  beam  theory  to  assume  that  a  certain 
portion  of  brackets  and  overlapping  webs  constitute  a 
rigid  region  on  the  ends  of  each  span  as  shown  in  Figure 
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3.   This  method  called  the  "Span  Point  Method"  has  been 
shown  by  finite  element  calculations  to  slightly  over 
extimate  the  extent  of  the  rigid  regions. 

To  calculate  the  length  of  the  rigid  section 
in  a  given  element  Lund  [6]  uses  the  expressions 

fj  =  S-p-j.  /  sin  a 

fJT  =   Sj   /   sin  a 
where  f   and  f    are  the  lengths  of  a  rigid  section  in 
beams  I  and  II  respectively  and  S   and  S^,.  are  the 
overlapping  web  depths  associated  with  each  beam.   The 
value  a  is  the  angle  of  intersection  between  the  two 
beams  as  shown  in  Figure  4a.   The  angle  is  measured  at 
the  vertex  of  the  overlapping  webs  on  the  beam  ends  and 
is  formed  by  the  lines  parallel  to  the  beam  axes  through 
that  vertex.   It  is  evident  that  in  the  most  common  case 
of  normal  beams  butting  flush  against  each  other,  the 
rigid  length  influence  on  each  beam  of  the  intersection 
is  the  depth  of  the  adjacent  beam.   Where  three  or  more 
beams  tie  in  at  the  same  point,  the  beams  with  the  greatest 
depth  govern  in  this  effect.   In  this  case  the  bracket 
influence  has  been  neglected. 

To  compute  the  influence  of  each  bracket  leg, 
the  approximation 

tj  =  d1   /Uh/d^  +  1} 


.b 
'1 

computed  rigid  length  in  the  bracket,  d-,  ,  is  the  bracket 


is  used  £ 6 J .   As  is  shown  in  figure  4b. ,  fc~  is  the 
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leg  length  and  h  is  the  beam  depth.   The  total  rigid 
length  of  one  end  of  a  given  beam  is  the  sum  of  the 
effects  of  the  bracket  and  the  beam  overlap  c- , 
Figure  4b. 

In  this  study,  when  bracket  influence  is 
considered,  intersecting  beams  will  be  normal  or  nearly 
normal  and  it  is  assumed  that  c,  =  fT. 

It  should  be  noted  here  that  the  preceeding 
expressions  are  based  on  known  values  of  geometry  which 
are  not  known  at  the  start  of  the  design  process. 
Iteration  from  an  assumed  starting  value  is  required. 

The  rigid  region  assumption  takes  on  special 
significance  in  the  limit  design  problem.   The  effect  is 
to  indirectly  specify  regions  within  the  structural  system 
in  which  plastic  hinges  may  not  form.   This  forces  span 
ends  toward  their  midpoints  while  the  total  external  work 
rate  part  of  the  virtual  work  relation  may  or  may  not 
remain  unchanged,  depending  upon  the  nature  of  the 
particular  mechanism  under  consideration.   In  any  event, 
the  formulation  of  the  virtual  work  relation  becomes 
significantly  more  complex  in  all  but  the  simplest  beam 
mechanism,  especially  in  the  external  work  part  of  the 
expression.   The  variability  of  rigid  length  values 
through  a  successive  iteration  steps  compounds  the 
problem.   Lund  [5]  side  steps  this  problem  by  calculating 
the  rigid  extensions  only  once  at  the  problem  beginning 
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and  holding  them  constant  thereafter. 

Rigid  section  examples  examined  in  this  study 
were  calculated  by  hand  with  rigid  extensions  assumed  to 
be  the  same  as  used  by  Lund. 
2.3   Combined  Mechanisms  and  Kinematic  Admissibility 

Two  final  considerations  are  very  important 
in  the  plastic  limit  design  calculations  performed  in 
this  study.  These  are  limitations  which  can  restrict 
combining  mechanisms  as  was  illustrated  in  section  2.1 
and  conflict  with  the  requirement  that  all  mechanisms 
be  kinematically  admissible. 

The  upper  bound  theorem  of  limit  design  requires 
that  all  velocity  fields  be  kinematically  admissible. 
This  is  a  requirement  that  deserves  special  attention  in 
formulating  the  limit  design  problem.   For  the  velocity 
field  in  a  plastic  collapse  mechanism  to  be  kinematically 
admissible  the  velocity  in  the  direction  normal  to  the 
axes  of  the  elements,  u  must  be  much  larger  than  in  line 
velocity,  w.   Figure  5  illustrates  that  only  slight 
increases  of  the  beam  lengths  are  required  to  move  point 
b  a  substantial  distance  laterally  through  beam  rotation 
to  point  b'  in  the  pictured  mechanism. 

Martin  [22]  states  that  the  hinges  of  joint 
rotation  mechanisms  must  form  at  the  joint  and  not  some 
finite  distance  out.   This  point  is  generally  overlooked 
in  the  literature,  but  is  of  utmost  importance  if  the 
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joint  mechanism  is  to  be  utilized  in  forming  combined 
mechanisms.   The  source  of  this  problem  is  the  requirement 
that  the  displacement  fields  of  all  limit  design 
mechanisms  be  kinematically  admissible.   If  a  multi- 
element joint  is  rotated  and  if  its  hinges  are  located 
anywhere  but  at  the  center  of  rotation,  kinematic 
admissibility  as  defined  previously  will  be  violated. 
Plastic  hinges  in  addition  to  the  joint  mechanism 
hinges  will  be  required  before  the  mechanism  can  rotate. 

The  assumption  of  a  rigid  bracket  region 
conflicts  with  the  use  of  joint  rotation  mechanisms 
since  the  rigid  regions  permit  hinges  at  their  extremes 
away  from  the  rotation  center.   Therefore,  while  the 
rigid  joints  can  certainly  be  used  as  components  of 
mechanisms,  they  cannot  be  considered  as  independent 
mechanisms.   This  all  but  eliminates  linear  combinations 
of  mechanisms  as  discussed  in  section  2.1,  if  joint 
rotations  are  required  to  form  the  combination.   When 
combined  mechanisms  are  formed  they  require  the 
formation  of  a  very  large  number  of  plastic  hinges 
with  correspondingly  high  values  of  internal  energy 
dissipation,  reducing  the  probability  that  such  a 
mechanism  would  contribute  to  collapse.   In  this  study, 
combined  mechanisms  were  considered  only  when  they 
consisted  of  seven  or  less  plastic  hinges. 

The  use  of  the  rigid  bracket  region  assumption 
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can  be  one  final  source  of  trouble.   If  the  rigid 
extensions  are  large  and  angles  between  hinge  points  are 
large,  the  proposed  mechanism  may  be  kinematically 
inadmissible.   In  Figure  6  which  approximates  the  rigid 
bracket  in  a  bilge  region,  u  and  w  are  of  the  same 
order  of  magnitude  and  the  mechanism  is  inadmissible 
as  is.   Were  hinge  c  moved  farther  to  the  right  the 
mechanism  could  become  admissible.   How  far  the  hinge 
must  be  moved  before  an  admissible  mechanism  will  be 
formed  is  a  question  of  judgement.   For  this  study,  it 
was  generally  required  that  the  height  of  the  elevated 
hinge,  h,  be  less  than  twice  the  distance  from  the 
corner  below  it  to  the  next  plastic  hinge  for  the 
mechanism  to  be  considered  kinematically  admissible. 
It  is  assumed  that  if  this  requirement  is  met,  the  value 
of  a  will  be  relatively  small. 

The  rigid  extension  assumption  is  likely  to  be 
of  greatest  importance  in  tanker  web  frames  which  do  not 
have  cross  ties.   These  frames  will  have  very  deep  webs 
and  correspondingly  large  overlap  regions.   One  frame 
of  this  type  was  examined.   In  frames  which  have  cross 
ties,  the  web  depths  will  be  less.   In  the  transverse 
frames  which  had  cross  ties,  the  rigid  extension 
assumption  was  not  used.   In  frames  of  this  type,  plastic 
hinges  were  permitted  anywhere  in  the  beam  element. 
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2.4   Summary 

In  this  study,  the  upper  bound  theory  of  limit 
design  is  used  to  develop  a  system  of  linear  relations 
which  must  be  satisfied  if  the  design  structure  is  to 
safely  carry  the  prescribed  loads.   A  linear  weighting 
or  object  function  is  formulated  based  on  the  lengths  of 
the  respective  structural  elements.   These  linear 
relations  combine  to  form  a  linear  programming  problem 
which  may  be  solved  to  find  the  minimum  required  plastic 
moments  of  the  various  structural  elements  in  the  system. 

It  is  assumed  that  the  plane  frameworks  of 
beams  are  made  of  materials  with  elastic,  perfectly 
plastic  moment  curvature  relations.   It  is  also  assumed 
that  transverse  shear  and  axial  forces  do  not  influence 
the  plastic  yielding  of  the  beams. 

Beam,  panel  and  fictitious  joint  rotation 
mechanisms  are  considered  in  the  formulation  of  linear, 
principal  of  virtual  work,  constraint  relations  as  well 
as  many  combined  mechanisms.   Combined  mechanisms  which 
are  very  complex  and  fundamental  sway  mechanisms  are 
not  considered. 

The  assumption  of  rigid  extensions  in  the 
area  of  beam  overlap  and  brackets  used  in  some  elastic 
designs  is  considered  in  one  case.   This  assumption 
disqualifies  the  use  of  fictious  joint  mechanisms  in 
combining  fundamental  mechanisms  on  grounds  of  kinematic 
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inadmissibility.   Mechanisms  involving  these  rigid 
regions  must  be  carefully  checked  for  kinematic 
admissibility.   The  rigid  extensions,  where  used,  are 
calculated  only  once  at  the  start  of  the  design 
sequence. 
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3.   OPTIMIZATION  TECHNIQUE 

3. 1   Linear  Programming 

It  was  stated  in  the  preceeding  chapter  that 
the  linear  equilibrium  constraint  relations  and  the 
minimum  weight  object  formed  a  linear  programming 
problem. 

The  general  linear  programming  problem  can 

be  described  as  follows 

Given  a  set  of  i  linear  inequalities  or 
equations  in  j  variables,  we  wish  to  find 
non-negative  values  of  these  variables 
which  will  satisfy  the  constraints  and 
maximize  or  minimize  some  linear  function 
of  the  variables  [30]. 

Mathematically  this  can  be  written  minimize 

(or  maximize)  the  linear  objective  function, 

N 
f  =  I    c,x  (3.1) 

j=l  3    1 

subject  to  the  linear  equality  constraints 

N 

I   a^ .x .  =  b.       i  =  1,....  (3.2) 


and 


x.  >  0  (3.3) 


The  values  of  a-.,  b.,  c.  are  assumed  to  be  known 
constants.   Other  notations,  expecially  vector  forms, 
also  appear  in  the  literature.   The  use  of  inequalities 
is  not  entirely  prohibited  in  equation  (3.2)  [29]. 
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the  use  of  slack  variables  permits  the  inequality 
N 


y   a.  .x,  > 

g    id  j  - 


b  (3.4) 

m 


to  be  expressed  as 
N 

I 

J 
XN+1  -  °  (3'6 


I    aijXj  '  XN+1  =  bm  <3-5> 


The  use  of  a  plus  sign  in  front  of  the  slack 
variable  would  signify  that  the  relation  was  of  the 
less  than  (<)  form. 

Any  set  of  x.  which  satisfies  the  constraints 
(3.2)  is  called  a  solution  to  the  linear  programming 
problem.   Any  solution  which  satisfies  the  non-negativity 
restrictions  of  equation  (3.3)  is  called  a  feasible 
solution.   Any  feasible  solution  which  optimizes  the 
objective  function  is  called  an  optimal  feasible  solution. 
Generally,  there  will  be  an  infinite  number  of  feasible 
solutions  to  a  linear  programming  problem.   From  all 
of  these  feasible  solutions,  one  must  be  found  which 
optimizes  the  objective  function. 

Linear  programming  problems  which  involve 
only  two  variables  can  be  solved  graphically.   The 
graphic  solution  to  these  two  variable  problems  can 
provide  a  great  deal  of  insight  into  the  more  general 
linear  programming  case  with  any  number  of  variables. 

In  Chapter  2,  a  simple  portal  frame  was 
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examined  and  a  system  of  linear  constraining  inequalities 

was  developed  which  gave  upper  bounds  on  the  plastic 

limit  moments  required  to  make  the  structure  safe. 

These  inequalities  in  two  variables  with  the  minimum 

weight  objective  function  form  a  linear  programming 

problem  which  is  illustrated  graphically  in  Figure  7. 

The  solid  lines  represent  the  various  constraint 

functions.   The  dotted  lines  are  the  parallel  lines  of 

constant  weight  of  the  object  function.   The  solid 

lines  with  cross-hatching  attached  bound  the  region  of 

feasible  solutions  (or  the  acceptable  design  region) . 

This  is  the  boundary  of  the  region  where  all  inequality 

constraint  relations  are  satisfied.   The  optimum  design 

is  found  at  the  point,  or  points,  where  the  object 

function 

Min  f  =  8M   +  4M, 
a     b 

has  its  lowest  value  while  just  touching  the  boundary  of 
feasible  solutions.  The  figure  shows  that  the  variables 
of  the  optimal  solution  are  approximately  14  =  1.8  and 

CI 

ML  ■  1.1.   To  find  the  exact  values,  simultaneous 

solution  of  the  equations  of  the  two  intersecting  lines, 

equations  (2.8a)  and  (2.10a)  is  required.   Restated, 

these  are 

2  M=  +  2  Mb  =  6 

11  M   +  3  M,  =  24 
a      Td 

Solution  of  these  equations  gives  exact  values 
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of  M  =  1.875  and  M,  =  1.125.   These  values  represent 
a  r> 

exact,  upper  bound  solutions  to  the  minimum  values  of 
plastic  limit  moment  required  to  make  the  design 
structure  of  Figure  1,  safe  under  the  prescribed  loads. 
Note  that  no  factors  of  safety  have  been  applied  at  this 
point. 

In  this  case  the  optimum  solution  was  defined 
to  be  at  single  vertex  formed  by  two  boundary  lines. 
It  is  also  possible  that  the  objective  function  could 
have  passed  through  two  vertices  as  well  as  all  points 
on  the  line  between  them.   In  cases  of  this  type,  it 
is  said  that  the  linear  programming  problem  has 
alternative  optima.   This  condition  occurs  when  the 
objective  one  of  the  extreme  points  is  optimal. 
Unfortunately  these  methods  do  not  tell  which  one  [29], 

Other  methods  must  be  used  to  solve  linear 
programming  problems.   There  is  presently  no  method 
available  which  will  find  an  optimal  solution  to  a 
linear  programming  problem  in  a  single  step.   All  of 
the  proceedures  for  solving  these  problems  are,  therefore, 
iterative.   If  the  number  of  constraints  is  large,  it 
becomes  very  difficult  to  find  any  feasible  solution 
and  still  more  difficult  to  find  an  optimal  solution. 
3.2   The  Revised  Simplex  Method 

The  revised  simplex  method  of  solving  linear 
programming  problems  is  a  revision  of  the  simplex  method 
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derived  by  Dantzig  in  1948.   The  term  "simplex"  has 
nothing  to  do  with  the  technique  as  it  is  now  used. 
To  properly  introduce  the  revised  simplex  method  it 
is  helpful  to  start  with  a  short  discussion  of  the 
simplex  method.   Both  techniques  are  algebraic 
iterative  procedures  which  will  solve  exactly  any 
linear  programming  problem  in  a  finite  number  of  steps, 
or  give  an  indication  that  the  solution  is  unbounded. 
There  are  a  finite  number  of  extreme  points  and  if 
there  is  an  optimal  solution,  one  of  the  extreme 
points  is  optimal.   The  simplex  method  is  a  technique 
which  moves  from  a  given  extreme  point  to  another 
adjacent  and  more  optimal  extreme  point,  along  the 
boundary  of  the  feasible  solution  region.   Of  all 
adjacent  extreme  points,  the  one  chosen  is  the  one 
that  gives  the  greatest  decrease  (or  increase)  in  the 
object  function.   At  each  iteration,  the  simplex 
method  defines  whether  or  not  the  associated  extreme 
point  is  optimal  and  if  not,  it  defines  the  next 
extreme  point  to  be  examined. 

Complete  discussions  on  the  simplex  and  revised 
simplex  methods  of  linear  program  are  by  nature  quite 
lengthy.   It  is  not  the  purpose  of  this  study  to 
examine  in  depth  the  mathematical  rational  for  either 
of  these  methods.   Detailed  information  on  theory, 
development  and  application  can  be  found  in  References 
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[29  -  31]. 

The  simplex  method  has  two  phases.   The  first 

phase  finds  a  first  basic  feasible  solution  which 

satisfies  the  condition  that  the  m  constraints  intersect 

and  that  there  are  only  m  nonzero  variables  x . .   These 

x.  variables  are  called  the  variables  of  the  basis. 
] 

The  remaining  (N  -  m)  of  the  variables  x .  are  equal 

] 

to  zero  and  called  non  basic.  At  the  end  of  the  first 
phase  the  results  are  structured  in  standard  canonical 
form.  In  order  to  satisfy  inequality  (3.3)  all  values 
of  b.  must  be  greater  than  zero  [8].  In  practice  this 
may  require  the  introduction  of  artificial  variables. 

When  the  canonical  form  of  the  first  phase  is 
obtained,  the  second  phase  begins.   Given  any  basic 
feasible  solution  an  iteration  step  is  made  toward  an 
optimal  feasible  solution  by  changing  a  single  vector 
in  the  basis.   In  the  revised  simplex  method,  the 
criteria  used  to  determine  the  vector  to  enter  and  leave 
the  basis  are  identical  with  the  simplex  method.   Changes 
in  the  basic  feasible  solution  always  lead  back  to  the 
canonical  form  with  the  variables  x  renumbered.   The 
essential  element  that  allows  progression  from  each 
basic  feasible  solution  is  the  explicit  knowledge  of  the 
representation  of  the  vectors  not  in  the  current  basis  in 
terms  of  the  basis  vectors. 

The  main  difference  between  the  simplex  method 
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and  the  revised  procedure  is  that  in  the  simplex 
method  all  elements  of  the  simplex  tableau  are  transformed 
by  means  of  elimination  formulas,  while  in  the  revised 
simplex  method  only  the  elements  of  an  inverse  matrix 
need  by  transformed  using  these  formulas. 

Since  its  development  the  revised  simplex 
method  and  in  particular  a  variation  which  employs  the 
product  form  on  the  inverse,  has  been  used  in  the  larger 
high-speed  computers  [29]. 

There  are  several  reasons  for  the  current 
widespread  use  of  the  simplex  method.   For  problems 
with  sparse  coefficient  matrices,  computation  is  reduced. 
The  revised  procedure  always  deals  with  original 
coefficients.   The  amount  of  new  information  required 
to  be  stored  is  less  than  in  the  simplex  method,  since 
only  the  inverse  and  the  solution  vector  must  be  stored 
and  not  the  complete  simplex  tableau  as  in  the  simplex 
method. 

To  formulate  the  computational  procedure  for 
the  revised  simplex  method  equations  (3.1)  through  (3.3) 
will  first  be  restated  in  vector  form.   In  the  general 
linear  programming  problem  a  minimum  value  of  the  object 
function 

c  X  (3.5) 

is  sought  subject  to 

A  X  =  b  (3.6) 
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and 

X  >  0  (3.7) 

letting  the  matrix  B  of  the  basis  of  m-dimensional  vectors 

correspond  to  the  first  m  vectors  of  A,  such  that 

B  XQ  =  b  (3.8) 

XQ  >  0  (3.9) 

where  Xq  =  (x-,fl,  x20'  ••••  X™  )  is  the  solution  vector. 

A  basic  feasible  solution  is  given  by 

X0  -  B"1  b  (3.10) 

The  linear  combinations  of  all  the  vectors  of 

A  in  terms  of  B  can  be  found  by 

X-j  =  ^-1  1-  j  =  1,2, .  .  .n  (3.11) 

where  X .  is  a  column  vector  and  P.  is  a  set  of  linearly 
-D  -1 

independent  vectors  associated  with  each  extreme  point 
and  can  be  expressed 

£j  -  xlj  El  +  x2j  ^2  +  ••••  +  Xmj  ^m     (3'12) 
The  initial  basic  feasible  solution  is  found 

by  setting  X.  =  P.. 
-0    -3 

For  a  given  iteration  solution  the  value  of 

the  objective  function  is 

z .  =  c,  x, .  +  c_  x_ .  +  ....  +  c   x  . 
3     1   lj     2   2j  m  mj 

j  =  1,  2,  ....,  n   (3.13) 

where  c.  are  cost  coefficients  of  those  vectors  in  the 

basis.   Converting  to  vector  form  z .  =  c„  X-  where  cn 

*  3-0-J-O 

=    (c,  ,  c-,  ....  c  )  is  a  row  vector. 
J.   a  m 

From  equation  (3.11)  z.  can  be  written  in  the 

3 
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form 

Z.  =  TT  P.  (3.14) 

J   — j 

where      7T  =  cQ  B"1  (3.15) 

tt_  is  an  m-dimensional  row  vector  called  the  multiplier 
vector.   The  individual  values  tt.  in  the  vector  are 
called  the  simplex  multipliers.   Given  the  value  of 
n  for  a  given  basis  B,  the  corresponding  z.  can  be 
calculated.   A  vector  P.  not  in  the  basis  is  evaluated 
by  computing  tt_  P  ■  -  c  ■  =  z  .  -  c  .  . 

The  information  needed  to  proceed  from  one 
feasible  solution  to  the  next  improved  feasible  solution 
is  the  original  data  consisting  of  A,  b  and  c  and 
equations  (3.10),  (3.11)  and  (3.15),  if  for  each  feasible 
basis  B,  B-i  is  known  explicitly  [29].   This  concept  is 
fundamental  to  the  revised  simplex  method  of  solving 
linear  programming  problems. 

The  revised  simplex  linear  programming  procedure 
has  been  applied  by  large  numbers  of  users  to  a  wide 
variety  of  problems  with  efficient  results.   Most  often, 
it  can  be  treated  as  a  "black  box"  operation  with  only 
general  knowledge  on  the  part  of  the  user.   It  forms  a 
part  of  most  computer  application  libraries  and  is 
available  at  little  or  no  increase  in  cost  over  the  basic 
computer  usage  fees  [8], 

The  subroutine  called  in  this  study  was  the 
Mathematical  Programming  System  Extended  (M.P.S.X.)  as 
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installed  on  the  IBM  370  Model  168  computer.   Specifics 
on  the  use  of  the  subroutine  are  outlined  in  the 
appendices. 
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4.   DISCUSSION  OF  RESULTS 

4  . 1   Single  Cross-Tie  Frames 

The  results  of  the  twenty-five  frames  evaluated 
in  Appendix  A  varied  significantly.   The  lightest  frame 
had  a  cross-tie  height  of  12.0  meters  (as  did  all  of 
Lund's  frames)  and  a  side  tank  breadth  of  12.0  meters. 
Lund's  lightest  frame  had  a  side  tank  width  of  16.0 
meters.   The  weights  of  the  lightest  limit  design  frame 
were  41.47  metric  tons  and  48.50  metric  tons  for  factors 
of  safety  of  1.5  and  2.0  respectively.   All  of  Lund's 
frames  are  based  upon  a  factor  of  safety  of  1.5  and  his 
lightest  structure  weighed  42.1  metric  tons. 

The  heaviest  limit  design  frame  structure  had 
a  cross-tie  height  of  10.0  meters  and  a  side  tank  width 
of  20.0  meters.   The  weights  of  this  frame  were  56.08 
metric  tons  and  63.45  metric  tons  for  factors  of  safety 
of  1.5  and  2.0  respectively. 

The  total  frame  weights  of  the  elastic  design 
structures  and  the  limit  design  structures  are  within 
10.0  metric  tons  for  all  frames  having  factors  of  safety 
of  1.5. 

The  distribution  of  frame  weights  resulting  from 
the  different  design  methods  is  significant. 

Member  proportions  were  kept  constant  through 
both  elastic  and  plastic  limit  design  optimization, 
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excluding  plate  flange  areas.   This  procedure  allows 
direct  comparison  of  the  relative  weights  of  members 
with  the  same  proportions  by  comparison  of  top  flange 
areas.   Comparisons  of  beam  proportions  are  made  using  a 
factor  of  safety  of  1.5  for  the  limit  design  results 
unless  stated  to  be  otherwise. 

The  limit  design  results  show  the  size  of 
Member  #5,  Figure  A.l,  to  be  larger  than  the  Lund  results 
for  almost  all  frames  considered.   The  limit  design 
results  with  factors  of  safety  of  1.5  gave  top  flange 
areas  of  124  square  centimeters  in  the  lightest  beam 
case  and  283  square  centimeters  in  the  heaviest  beam 
case.   Lund's  heaviest  lower  side  tank  beam  had  a  top 
flange  area  of  127  square  centimeters  and  his  lightest 
had  an  area  of  65  square  centimeters.   As  might  be 
expected  the  heaviest  beams  were  those  with  the  longest 
spans.   The  limit  design  lower  side  tank  beams  have 
roughly  twice  the  cross  sectional  area  of  the  elastic 
designs  with  the  same  side  tank  breadth  in  all  cases. 
For  example,  Lund's  frames  with  16  meter  side  tank 
widths  gave  values  of  92  and  127  square  centimeters  for 
top  flange  area.   For  similar  frames  the  limit  design 
results  gave  values  of  198,  200,  203,  205,  and  207  for 
top  flange  areas  of  the  lower  side  tank  beams.   The 
dimensions  of  the  Lund  elastic  frames  may  be  found  in  the 
second  rows  of  Tables  3  and  4  of  Reference  [5 J . 
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The  model  used  in  the  limit  design  study 
considered  all  loads  acting  on  the  structure.   Lund 
does  not  discuss  his  model  is  detail.   It  is  not  known 
if  he  accounted  for  the  influence  of  loads  on  the  rigid 
extensions  in  his  beam  analysis.   If  loads  in  the  rigid 
regions  were  neglected  in  the  elastic  designs,  this 
could  account  for  his  smaller  required  size  for  the 
lower  side  tank  beam. 

With  one  notable  exception  other  beam  elements 
were  generally  slightly  lighter  for  the  limit  design 
case.   The  values  calculated  for  the  dimensions  of 
Element  #7  were  much  less  in  the  limit  design  results. 
Representative  values  of  from  69  to  101  square  centimeters 
were  calculated  for  the  top  flange  areas  of  the  limit 
design  frames  versus  a  range  of  134  to  253  for  the 
elastic  designs.   In  general,  the  elastic  designs  were 
roughly  twice  the  weight  of  the  plastic  limit  designs 
for  this  beam  element  (i.e.  Element  #7).   This  difference 
may  have  again  been  caused  by  unknown  differences  in 
models. 

All  comparisons,  so  far,  have  been  based  on 
identical  values  of  factors  of  safety  of  1.5  for  elastic 
and  plastic  limit  designs.   Figure  A. 5  shows  that  when 
a  factor  of  safety  of  2.0  is  applied  to  the  plastic  limit 
designs,  the  total  frame  weights  become  higher  than  all 
but  one  of  the  frames  designed  by  elastic  methods.   The 
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effect  of  increases  in  factors  of  safety  on  individual 
beam  elements  is  illustrated  in  Figure  A. 4  for  one  beam 
element  for  the  various  plastic  limit  design  frame 
geometries.   Values  of  net  frame  area  A   are  plotted 
against  factors  of  safety  for  Member  #1  of  all  frames 
evaluated.   The  dots  represent  values  for  the  elastic 
designs  of  Lund.   The  increase  in  net  area  is  seen  to 
be  nearly  linear.   On  this  figure  only  the  line  representing 
frames  having  cross-tie  heights  of  12.0  meters  should  be 
compared  with  the  elastic  design  values  since  all  elastic 
designs  had  12.0  meter  cross-tie  heights. 

Plots  for  other  beam  elements  would  be  similar 
to  Figure  A. 4,  but  there  could  be  up  to  25  different 
lines  if  each  of  the  frames  considered  had  different 
values  of  AT  for  the  element  considered. 

As  would  be  expected  from  looking  at  Figure 
A. 5,  the  weights  of  the  individual  limit  design  beam 
elements  are  in  general  heavier  than  Lund's  results 
when  a  factor  of  safety  of  2.0  is  applied  to  the  limit 
design  results. 

The  design  effort  required  in  evaluating  the 
twenty-five  frames  discussed  here  was  minimal.   Total 
computation  cost  was  $91.65.   This  results  in  a  computation 
cost  per  frame  of  $3.66.   The  limit  design  problem  was 
greatly  simplified  by  using  the  assumption  that  the 
influence  of  brackets  and  beam  overlap  could  be  neglected. 
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Lund  did  not  use  this  assumption.   His  elastic  designs 
all  assumed  rigid  joint  regions.   The  elastic  design 
computation  time  is  much  larger  than  for  the  plastic 
designs.   Lund  [5]  cites  CPU-time  of  50  to  60  seconds 
on  a  UNIVAC  1108  for  design  of  one  frame  structure.   The 
limit  design  solutions  of  Appendix  A  required 
approximately  8  seconds  in  CPU-time  for  one  frame 
structure  on  the  IBM  370  including  M.P.S.X.  linear 
programming  solutions. 
4 . 2   Zero  Cross-Tie  Frame  with  Rigid  Corner  Assumption 

The  results  obtained  in  Appendix  B  were  for  a 
single  frame  based  on  the  geometry  of  Lund's  frame 
numbered  OC22,  49. 2t.   The  values  of  total  frame  weight 
for  a  number  of  factors  of  safety  are  plotted  in  Figure 
A. 6.   The  weight  of  the  elastic  parent  frame  is  plotted 
as  a  point  at  a  factor  of  safety  of  1.5.   From  the  figure 
it  may  be  seen  that  the  weight  of  the  limit  design  frame 
calculated  with  a  factor  of  safety  of  1.5  is  twenty  per 
cent  lighter  than  the  elastic  frame.   At  a  factor  of  safety 
of  2.0,  the  limit  design  structure  weighs  approximately 
one  half  ton  less  than  the  elastic  design  structure. 

In  the  top  line  of  Table  4,  Reference  [5]  ,  Lund 
lists  the  dimensions  of  his  elastic  design  frame.   Table 
A.  2. 4  lists  values  for  the  limit  design  frame  for  factors 
of  safety  of  1.0  to  2.5  in  0.5  increments.   If  the  limit 
design  dimensions  with  a  factor  of  safety  of  1.5  are 
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compared  to  the  Lund  results,  it  will  be  seen  that  the 
lower  side  tank  girder,  Member  #2,  is  again  substantially 
larger  in  the  limit  design  results  than  in  the  elastic 
design  results.   Using  top  flange  areas  as  a  measure  of 
relative  weight,  the  limit  design  beam  has  a  value  of 
192  cm  ,  while  the  elastic  design  beam  has  an  area  of 
167  cm  .   This  gives  a  relative  difference  in  beam  weights 
of  approximately  15  per  cent. 

As  mentioned  previously,  it  is  not  known  if  Lund 
accounted  for  the  influence  of  loads  on  the  rigid 
extensions.   The  limit  design  calculations  accounted  for 
all  loads  on  the  structure  including  the  rigid  regions. 
If  Lund  did  not  do  this  in  his  design  model,  his  results 
could  be  lighter  than  the  plastic  limit  design  values. 

With  the  exception  of  Member  #5,  in  the  side  tank 
top,  all  of  the  remaining  limit  design  beam  elements  were 
about  half  the  weight  of  their  elastic  design  counterparts. 
Member  #5  in  the  side  tank  top  was  about  one  fourth  the 
weight  of  its  elastic  design  counterpart. 

The  limit  design  calculations  done  in  Appendix 
B  were  done  by  hand  as  opposed  to  Lund's  computer  approach. 
The  hand  operations  were  performed  to  allow  close  inspection 
of  each  mechanism.   The  results  were  taken  to  the  computer 
for  solution  of  the  linear  programming  problem. 
Computation  time  was  approximately  6  seconds.   Again,  the 
CPU-time  for  Lund's  solutions  averaged  50-60  seconds  per 
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frame  structure.   The  time  required  to  do  the  limit 
design  hand  calculations  was  approximately  30  hours  on 
the  first  time  through.   This  time  could  be  reduced  to 
about  four  hours  for  successive  runs  by  using  the  models 
already  developed. 

These  solutions  could  also  be  adapted  for 
computer  solution  as  was  done  in  Appendix  A  using  a  model 
similar  to  the  one  used  in  the  hand  calculations  and  a 
heavy  reliance  upon  computer  conditional  tests.   The 
conditional  tests  would  be  required  to  check  the  possible 
mechanisms  for  kinematic  admissibility. 
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5.   SUMMARY  AND  CONCLUSIONS 

An  analytical  study  has  been  undertaken  to 
investigate  the  use  of  the  upper  bound  theorem  of  limit 
design  in  the  design  of  transverse  web  frames  in  ships. 
Frames  have  been  developed  using  limit  design  with 
results  compared  to  frames  of  elastic  design. 

In  Appendix  A  the  influence  of  beam  overlap  and 
brackets  was  neglected.   This  assumption  was  not  made  in 
the  elastic  design  case.   The  results  based  on  this 
assumption  were  frame  structures  of  roughly  equal 
weight  to  the  elastic  frame  designs  with  some  frames 
less  than  the  least  weight  elastic  designs  if  the  same 
factor  of  safety  is  applied  in  both  methods.   For  a  factor 
of  safety  2.0  the  plastic  limit  design  structures  are 
approximately  10  per  cent  heavier. 

It  is  significant  that  while  the  overall 
structural  weights  are  similar  in  this  comparison,  the 
distribution  of  weight  is  somewhat  different.   Frame 
structures  designed  by  elastic  methods  have  their 
largest  members  in  one  location  while  structures  designed 
by  the  plastic  method  have  their  largest  members  in 
another  location. 

In  Appendix  B  the  effect  of  beam  overlap  and 
bracket  influence  was  idealized  by  assuming  portions  of 
these  regions  to  be  infinitely  rigid. 
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Structural  weight  calculated  using  this  assumption 
was  approximately  equal  to  the  elastic  design  if  a  factor 
of  safety  of  2.0  was  applied  to  the  limit  design  solution. 
If  a  factor  of  safety  of  1.5  is  applied  to  the  limit  design 
solution,  its  result  gives  a  twenty  per  cent  reduction  of 
weight  over  the  elastic  design. 

A  specific  factor  of  safety  was  not  specified 
for  the  plastic  designs  but  rather  a  range  of  values  was 
examined.   The  important  advantage  of  limit  design  is  that 
it  allows  the  designer  to  define  exactly  the  collapse 
moments  for  a  given  set  of  loading  conditions.   Factors 
of  safety  may  be  applied  to  these  known  quantities  as 
required. 

In  the  case  of  the  elastic  design  methods, 
factors  of  safety  are  defined  somewhat  arbitrarily 
without  specific  knowledge  of  when  collapse  will  occur. 
The  knowledge  of  specific  collapse  moments  allows 
structures  designed  using  plastic  design  methods  to  be 
much  more  efficient. 

The  geometry  of  individual  beam  cross  sections 
was  not  optimized  or  allowed  to  vary  from  the  Lund 
proportions  in  this  study.   This  was  to  permit  direct 
comparison  with  the  Lund  designs.   The  limit  design 
solutions  could  possibly  have  been  made  much  lighter  by 
use  of  more  efficient  beam  cross  section  geometry. 

In  the  case  where  the  assumption  of  rigid 
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corners  was  made,  no  attempt  was  made  to  calculate 
structural  features  within  the  rigid  assumption  zone. 
If  this  assumption  is  to  be  used,  corner  scantlings  will 
have  to  be  developed  by  separate  analysis. 

Neglecting  other  possible  advantages,  the 
design  costs  using  plastic  limit  design  are  very  cheap. 
A  design  solution  for  a  rather  sophisticated  tanker  web 
frame  can  be  had  for  very  few  dollars.   This  allows  a 
number  of  different  geometries  to  be  examined  as  in 
Appendix  A  for  the  price  of  a  single  design  found  by 
elastic  methods. 
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Figure  1.   One  bay  portal  frame. 
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Figure  2.   Fundamental  collapse  mechanisms 
of  portal  frame  in  Figure  1. 
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Figure  3.  Tanker  transverse  web  frame  model 
with  rigid  extensions  due  to  beam 
overlap  and  bracket  influence. 


-50- 


Figure  4a. 


Corner  model  for  calculation 
of  rigid  extensions  due  to  beam 
overlap. 


Figure  4b.   Corner  model  with  bracket  for 

calculation  of  rigid  extensions. 
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Figure  5.   Kinematically  admissible  collapse 
mechanism. 


Figure  6.   Kinematically  inadmissible 
collapse  mechanism  with  w 
approximately  equal  to  u. 
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Figure  7.   Graphical  solution  to  linear 
programming  problem  of  portal 
frame  example. 
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Figure  A.l.   Tanker  web  frame  model  with 
one  cross-tie. 
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Figure  A. 2.   Loading  conditions  considered  in  tanker 
web  frame  design. 
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Figure  A. 3.   Tanker  web  frame  model  without 
cross-ties  and  with  rigid 
extensions. 
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Figure  A. 4 


Frame  cross  sectional  area  vs.  factor  of 
safety  for  beam  Element  #1  of  tanker  web 
frame  with  a  single  cross-tie. 
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Figure  A. 5 


Various  frame  weights  for  limit  design  and 
Lund  elastic  frames  having  a  single  cross- 
tie  located  twelve  feet  above  the  base  line, 
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Figure  A. 6 


Weights  of  zero  cross-tie  tanker  web  frame 
with  rigid  joints  for  various  factors  of 
safety. 
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APPENDIX  A  -  TANKER  FRAME  COMPARISON 

Modern  supertankers  are  designed  with  large 
variations  from  ship  to  ship  in  the  internal  primary 
stiffening  system.   The  number  of  cross-ties  varies  from 
zero  to  three.   Structures  may  be  designed  with  or  without 
longitudinal  side  girders.   Bottom  center  girders  may 
be  low  or  high.   When  compared  on  the  basis  of  weight 
to  strength  ratios,  some  designs  should  be  superior  to 
others  [5j. 

The  present  study  comprises  investigations 
of  variations  of  a  single  transverse  web  frame 
configuration . 

A  FORTRAN  computer  program  was  developed  to 
evaluate  the  limit  moment  values  of  a  tanker  webframe 
having  a  single  cross-tie  located  in  the  wing  tank. 
The  vessel  considered  is  identical  to  the  one  considered 
by  Lund  in  Reference  [5],   The  vessel  is  a  typical 
supertanker  with  the  following  principal  dimensions: 
Length  L^   =  320.0  m 

Breadth  B  =   5  0.0  m 

Depth  D  =   25.0  m 

Draft  d  =   19.5  m 

Camber  (all  in  side  tank)  =  1.0  m 
Spacing  of  transverse  frames  =  5.0  m 
The  basic  geometry  of  the  frame  studied  is 
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shown  in  Figure  A.l.   The  program  was  developed 
specifically  to  evaluate  variations  of  this  frame,  but 
it  could  be  made  more  general  without  extensive 
modification.   As  written,  the  required  input  values 
are  the  cross-tie  height  above  the  base  line  S2  and 
the  wing  tank  breadth  S5,  both  values  being  in 
centimeters.   All  other  values  are  either  written 
into  or  calculated  by  the  program  as  needed,  based 
on  the  parent  vessel  geometry  and  the  four  loading 
conditions  prescribed  by  Lund  [5]  and  to  be  described 
later  in  this  appendix.   System  lines  are  assumed  to 
be  formed  by  the  inside  surface  of  the  shell  plating 
and  the  center  tank  face  plating  of  the  longitudinal 
bulkhead.   It  was  found  when  beams  were  sized  that 
plastic  neutral  axes  for  the  various  beam  elements 
remained  at  or  very  near  these  system  lines. 

The  possible  influence  of  beam  overlap  and 
brackets  was  neglected  with  plastic  hinges  allowed  to 
form  anywhere  along  the  lengths  of  individual  beam 
elements.   Combined  plastic  collapse  mechanisms  which  are 
very  complex  and  are  unlikely  to  be  actual  collapse 
mechanisms  have  been  eliminated  from  computation 
a  priori  to  reduce  computing  costs. 

The  range  of  cross-tie  heights  examined  is 
from  10  meters  to  14  meters  in  one  meter  intervals. 
The  range  of  side  tank  breadths  examined  is  from  12 
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meters  to  20  meters  in  2  meter  intervals.   In  all  a  total 
of  twenty-five  separate  variations  of  the  same  single 
cross-tie  web  frame  were  examined. 

The  loading  conditions  applied  were  identical 
to  those  used  by  Lund  in  his  elastic  designs.   The 
various  conditions  are  shown  in  Figure  A. 2.   The  loading 
models  as  applied  represent  the  following  situations 
[5J. 

1.  Empty  center  tank,  cargo  in  side  tanks,  full  draft 
with  dynamic  additions. 

2.  Empty  side  tanks,  cargo  in  center  tank,  full  draft 
with  dynamic  additions. 

3.  Empty  side  tanks,  ballast  in  center  tank,  draft 
of  0.25  D. 

4.  Ballast  in  center  and  side  tanks,  draft  of  0.35  D. 

As  stated  previously,  the  program  input 
variables  are  cross-tie  height  and  side  tank  breadth, 
S2  and  S5  respectively.   The  input  format  is  free 
with  the  real  values  of  S2  and  S5  in  centimeters  listed 
in  the  same  order  shown  here.   The  program  sequence 
exits  when  a  zero  value  for  S2  is  read  so  the  last 
data  card  should  have  real  values  of  zero  and  any  other 
real  constant  listed  consecutively. 

The  output  values  of  the  FORTRAN  program 
include  the  coefficients  of  the  internal  energy 
dissipation  portion  of  the  limit  design  upper  bound 


-62- 


theorem.   These  coefficients  are  presented  as  the 
values  a. .  of  the  coefficient  matrix  A  in  the  limit 
design  problem  constraint  relation.   The  associated 
values  of  the  external  work  rates  in  the  principle  of 
virtual  work  relations  are  also  given  as  output.   These 
work  values  represent  the  right  hand  side  of  the  constraint 
equations  b.  . 

The  output  information  is  all  punched  on 
cards  to  be  used  as  input  data  in  the  linear  programming 
problem  to  be  solved  by  the  M.P.S.X.  subroutine. 
The  format  of  the  punched  cards  is  quite  important, 
but  is  not  discussed  here.   A  discussion  of  M.P.S.X. 
input  format  is  given  in  the  numerical  routine  appendix. 

Finally,  a  punched  card  is  inserted  at  the 
front  of  each  output  data  deck  listing  the  values  of 
S2  and  S5.   This  is  to  enable  the  user  to  identify 
the  various  output  data  sets  when  more  than  one  frame 
is  examined. 

The  program  was  checked  for  logical  and 
numerical  errors  by  comparing  output  values  and  hand 
calculation  results  on  the  web  frame  having  S2  =  12.0 
meters  and  S5  =  12.0  meters. 

The  complete  listing  of  the  program  with 
sample  data  values  is  presented  on  the  pages  that 
follow. 
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The  output  data  from  the  program  is  not  listed 
in  this  document  as  computed.   To  list  this  information 
would  require  a  separate  volume. 

To  complete  the  M.P.S.X.  linear  programming 
data  set,  objective  function  cards  listing  the  various 
beam  lengths  are  added  to  the  FORTRAN  program  output 
cards.   Finally  the  desired  M.P.S.X.  control  language 
is  inserted  and  the  problem  run  to  get  optimum  upper 
bound  values  for  the  plastic  limit  moment  for  each 
beam  in  each  frame.   It  should  again  be  noted  that  the 
very  rigid  M.P.S.X.  input  format  requirements  outlined 
in  the  numerical  routine  appendix  must  be  followed 
precisely.   Failure  to  do  so  may  result  in  the  linear 
programming  routine  being  abandoned  or  worse,  in 
erroneous  results. 

In  order  to  compare  the  web  frame  sizes  required 

by  the  plastic  limit  design  problem  to  Lund's  elastic 

design  results,  the  values  of  plastic  limit  moment  are 

used  to  size  beam  elements  based  on  the  beam  proportions 

used  by  Lund  in  his  analysis.   A  0.2  cm  corrosion 

allowance  is  then  applied  to  plating  thicknesses.   Lund 

based  his  beam  proportions  on  a  set  of  starting  values 

which  he  scaled  up  to  meet  stress  control  requirements. 

The  values  of  plate  flange  area  A   are  likely  to  be 

P 

drived  by  design  factors  not  considered  here  and  are 
assumed  to  be  held  constant.   Lund  does  not  list  thickness 
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values  for  plate  flanges.   It  is  therefore  assumed  that 
the  full  area  of  A  may  be  considered  to  be  part  of 
the  beam  to  be  sized.   The  starting  values  of  web 
depth,  web  thickness,  top  flange  area  and  plate  flange 
area  as  used  by  Lund  and  in  this  study  for  scaling  beams 
is  listed  in  Table  A. 1.1. 

Lund  assumed  his  designs  to  be  of  mild  steel. 

2 
Stress  values  were  restricted  to  1500  kp/cm  versus  the 

o 
yield  stress  value  of  2249.8  kp/cm  .   This  gives  an 

effective  factor  of  safety  of  1.50.   There  is  some 

question  as  to  the  "proper"  factor  of  safety  to  be 

applied  to  the  values  of  plastic  limit  moment  to  make 

a  reasonable  comparison  between  the  elastic  and  plastic 

designs. 

Values  of  beam  geometry  for  the  plastic  limit 
designs  were  calculated  using  factors  of  safety  of  1.5, 
2.0,  and  2.5.   The  dimensions  calculated  included  those 
listed  in  the  heading  of  Table  A. 1.1,  as  well  as  the 
total  cross  sectional  area  of  the  beam  excluding  plate 
flange  area.   This  total  cross  sectional  area  will 
henceforth  be  called  A_. 

In  calculating  the  values  named  above  the 
following  relationships  were  satisfied. 


where 


M„  =  a     z 

0   y  p 


zp  =  A  a 
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In  these  relationships  M0  is  the  plastic  limit 
moment  calculated  in  the  limit  design  problem,  o   is 
the  material  yield  stress,  A  is  the  total  required 
area  of  the  cross  section  and  a  is  the  distance  between 
centroids  of  the  two  equal  areas  lying  above  and  below 
the  plastic  neutral  axis. 

The  values  A™,  web  depth  h,  web  thickness  t 
and  top  flange  area  A„,,  are  listed  in  the  tables  at  the 
end  of  this  appendix,  starting  with  Table  A. 1.2.   The 
numbers  in  the  frame  column  represent  values  of  S2 
and  S5  respectively.   For  example  "12  -  14"  denotes  that 
this  is  the  frame  where  S2  =  12.0  meters  and  S5  =  14.0 
meters.   Dimensions  with  factors  of  safety  of  1.5  and 
2.0  are  listed  for  all  beam  elements  exclusive  of 
Element  #4.   It  was  found  that  the  limiting  mechanisms 
of  the  limit  design  problem  never  included  this  element. 
Those  mechanisms  in  which  failure  of  the  member  might 
have  occurred  had  high  values  of  internal  energy 
dissipation  with  low  external  work  rates  when  compared 
with  the  other  possible  collapse  mechanisms. 

It  is  known  that  Member  #4  is  in  compression 
and  subject  to  buckling.   Plastic  buckling  is  not 
considered  in  this  study.   In  order  to  provide  some 
measure  of  comparison,  the  dimensions  found  by  Lund 
in  his  study  are  used  unaltered  in  this  analysis  when 
final  frame  weights  are  computed.   In  these  computations 
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Element  #4  is  assumed  to  have  a  cross  sectional  area  of 

o 
432  cm  with  a  web  height  of  153  cm  for  all  cases. 

Figure  A. 4  illustrates  variations  of  net 
frame  area  A   for  various  factors  of  safety  as  calculated 
for  Member  #1. 

Total  web  frame  weights  were  calculated  for 

one  half  of  the  complete  frames.   To  be  consistent  with 

Lund,  comparisons  are  based  on  the  assumption  that  the 

relative  amount  of  local  stiffening  and  relative  bracket 

sizes  are  equal  for  all  investigated  frames  and  these 

components  are  neglected.   Overlaps  are  compensated 

for  by  subtracting  one  half  of  the  average  weight  in 

the  joint  region  for  each  beam.   Total  web  frame  weights 

were  therefore  calculated  using  the  following  relation: 

9 
W  =  Pj^H  AT  -  h    (C.  +  Ck)  AT} 

where  p  is  material  density,  1-  is  the  member  length 
and  c.  and  c,  are  the  overlaps  at  the  respective  beam 

J        K 

ends.   This  relationship  was  also  used  by  Lund  [5]. 

The  weight  in  metric  tons  of  each  frame 
considered  is  listed  in  Table  A. 1.18.   Weights  were 
calculated  based  on  factors  of  safety  of  1.5  and  2.0. 

Lund  assumed  a  constant  cross-tie  height  of 
12.0  meters.   The  limit  design  frame  group  with  the 
12.0  meter  cross-tie  height  had  the  lowest  weights, 
as  a  group.   Figure  A. 5  shows  values  of  frame  weight 
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versus  side  tank  breadth  for  single  cross-tie  frames 
with  cross-tie  heights  of  12.0  meters.   Solid  lines 
represent  limit  design  results.   Dashed  lines  represent 
all  values  for  single  cross-tie  frames  given  by  Lund. 
There  are  two  curves  joining  the  elastic  frame  weights. 
The  one  to  the  upper  right  represents  weight  values 
where  the  bulkhead  girders  are  outboard  of  the 
longitudinal  bulkhead.   The  remaining  elastic  design 
curve  gives  weight  values  for  frames  with  bulkhead 
girders  inboard  of  the  longitudinal  bulkhead.   The  two 
limit  design  curves  both  assume  the  bulkhead  girders  to 
to  inboard  of  the  longitudinal  bulkhead.   The  two 
curves  give  weights  for  factors  of  safety  of  1.5  and 
2.0  for  the  various  frames. 
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TABLE  A. 1.1 

Initial  Member  Dimensions  on  which 
Beam  Proportions  are  Based 


Member 
No. 

Height 
(cm) 

Web 
Thickness 
(cm) 

Top  Flange 
Area 
(cm2) 

Plate  Flange 
Area 
(cm2) 

1 

300 

1.5 

120 

660 

2 

300 

1.5 

120 

660 

3 

300 

1.4 

80 

750 

4 

160 

2.0 

150 

150 

5 

400 

1.8 

120 

750 

6 

300 

1.5 

120 

420 

7 

300 

1.5 

120 

510 

8 

240 

1.25 

80 

750 

9 

400 

2.0 

120 

750 
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Table  A. 1.2  -  Dimensions  of  Member  #1 
(Factor  Safety  =  1.5) 


Frame 

AT, 

h 

t 

AF 

(S2    -   S5) 

(cm' 

!) 

(cm) 

(cm) 

(cm2) 

10-12 

428. 

93 

244. 

63 

1, 

.42 

79.79 

10-14 

ii 

n 

n 

ii 

10-16 

i) 

ii 

ii 

ii 

10-18 

ii 

ii 

M 

M 

10-20 

ii 

ii 

H 

ii 

11-12 

388. 

41 

232. 

04 

1 

.36 

71.79 

11-14 

ii 

ii 

ii 

ii 

11-16 

M 

n 

ii 

ii 

11-18 

n 

ii 

H 

ii 

11-20 

ii 

M 

ii 

ii 

12-12 

348. 

82 

219. 

09 

1 

.30 

64.00 

12-14 

n 

n 

H 

ii 

12-16 

ii 

n 

ii 

ii 

12-18 

■I 

ii 

M 

n 

12-20 

ii 

ii 

ii 

ii 

13-12 

311. 

22 

206. 

09 

1, 

.33 

56.63 

13-14 

ii 

ii 

M 

n 

13-16 

■1 

ii 

ii 

n 

13-18 

ii 

M 

ii 

M 

13-20 

ii 

n 

n 

n 

14-12 

309. 

06 

205. 

32 

1 

.33 

56.21 

14-14 

305. 

83 

204. 

17 

1 

.32 

55.58 

14-16 

275. 

61 

193. 

04 

1 

.17 

49.68 

14-18 

n 

ii 

n 

n 

14-20 

ii 

ii 

ii 

ii 

-79- 


Table  A. 1.3  -  Dimensions  of  Member  #2 
(Factor  Safety  =  1.5) 


Frame 

AT0 

h 

t 

AF 

(cm2) 

(S2  -  S5) 

(cm2) 

(cm) 

(cm) 

10-12 

613.10 

295.49 

1.68 

116.42 

10-14 

642.58 

302.88 

1.71 

122.32 

10-16 

669.92 

309.58 

1.75 

127.79 

10-18 

699.33 

316.64 

1.78 

133.68 

10-20 

728.71 

323.55 

1.82 

139.58 

11-12 

658.35 

306.77 

1.73 

125.47 

11-14 

687.78 

313.89 

1.77 

131.37 

11-16 

718.22 

321.10 

1.81 

137.47 

11-18 

748.63 

328.15 

1.84 

143.58 

11-20 

782.15 

335.76 

1.88 

150.32 

12-12 

701.43 

317.14 

1.79 

134.11 

12-14 

733.95 

324.77 

1.82 

140.63 

12-16 

766.44 

332.22 

1.86 

147.16 

12-18 

799.95 

339.74 

1.90 

153.89 

12-20 

833.42 

347.09 

1.94 

160.63 

13-12 

745.49 

327.43 

1.84 

142.95 

13-14 

780.06 

335.29 

1.88 

149.89 

13-16 

813.55 

342.74 

1.91 

156.63 

13-18 

848.05 

350.26 

1.95 

163.58 

13-20 

884.61 

358.06 

1.99 

170.95 

14-12 

771.68 

333.40 

1.87 

148.21 

14-14 

811.46 

342.28 

1.91 

156.21 

14-16 

859.55 

352.73 

1.96 

165.89 

14-18 

897.14 

360.70 

2.00 

173.47 

14-20 

932.61 

368.07 

2.04 

180.63 

-80- 


Table  A. 1.4  -  Dimensions  of  Member  #3 
(Factor  Safety  =  1.5) 


Frame 

(cm2) 

h 

t 

(cm2) 

(S2   -   S5) 

(cm) 

(cm) 

10-12 

426.89 

259.46 

1.41 

59.84 

10-14 

ii 

n 

ii 

n 

10-16 

■I 

ii 

ii 

ii 

10-18 

n 

it 

■■ 

ii 

10-20 

ii 

n 

ii 

n 

11-12 

386.19 

245.93 

1.35 

53.76 

11-14 

ii 

ii 

ii 

ii 

11-16 

ii 

ii 

ii 

■i 

11-18 

■i 

H 

ii 

ii 

11-20 

n 

H 

ii 

H 

12-12 

346.40 

231.99 

1.28 

47.84 

12-14 

ii 

n 

ii 

n 

12-16 

ii 

H 

ii 

ii 

12-18 

■I 

ii 

ii 

ii 

12-20 

ii 

ii 

ii 

n 

13-12 

309.68 

218.40 

1.22 

42.40 

13-14 

ii 

ii 

ii 

ii 

13-16 

ii 

ii 

ii 

H 

13-18 

ii 

n 

ii 

ii 

13-20 

ii 

ii 

ii 

ii 

14-12 

113.03 

125.14 

.78 

13.92 

14-14 

135.76 

138.78 

.85 

17.12 

14-16 

273.87 

204.35 

2.15 

37.12 

14-18 

ii 

ii 

ii 

ii 

14-20 

ii 

ii 

ii 

ii 

-81- 


Table  A. 1.5  -  Dimensions  of  Member  #5 
(Factor  Safety  =  1.5) 


Frame 

(cm2) 

h 

t 

AFo 

(S2  -  S5) 

(cm) 

(cm) 

(cm2) 

10-12 

956.60 

408.01 

2.04 

124.86 

10-14 

1216.50 

462.50 

2.28 

160.43 

10-16 

1488.72 

513.62 

2.51 

197.86 

10-18 

1766.22 

561.12 

2.73 

236.14 

10-20 

2048.14 

605.69 

2.93 

275.14 

11-12 

980.67 

413.35 

2.06 

128.14 

11-14 

1237.32 

466.60 

2.30 

163.29 

11-16 

1508.43 

517.13 

2.53 

200.57 

11-18 

1783.80 

564.00 

2.74 

238.57 

11-20 

2064.64 

608.20 

2.94 

277.43 

12-12 

1003.68 

418.39 

2.08 

131.29 

12-14 

1257.09 

470.46 

2.32 

166.00 

12-16 

1526.05 

520.26 

2.54 

203.00 

12-18 

1800.34 

566.69 

2.75 

240.86 

12-20 

2079.08 

610.39 

2.95 

279.43 

13-12 

1022.49 

422.46 

2.10 

133.86 

13-14 

1274.78 

473.89 

2.33 

168.43 

13-16 

1541.60 

523.00 

2.55 

205.14 

13-18 

1813.78 

568.88 

2.76 

242.71 

13-20 

2092.48 

612.41 

2.96 

281.29 

14-12 

1099.77 

438.83 

2.17 

144.43 

14-14 

1338.20 

485.99 

2.39 

177.14 

14-16 

1556.11 

525.54 

2.56 

207.14 

14-18 

1827.21 

571.05 

2.77 

244.57 

14-20 

2104.85 

614.27 

2.96 

283.00 

-82- 


Table  A. 1.6  -  Dimensions  of  Member  #6 
(Factor  Safety  =1.5) 


Frame 

AT^ 

h 

t 

AF. 

(S2    -   S5) 

(cmz 

!) 

(cm) 

(cm) 

(cm2) 

10-12 

384. 

13 

230.67 

1.35 

70.95 

10-14 

■I 

ii 

ii 

n 

10-16 

H 

•■ 

M 

ii 

10-18 

n 

ii 

ii 

ii 

10-20 

ii 

M 

M 

■I 

11-12 

342. 

38 

216.92 

1.28 

62.74 

11-14 

■i 

ii 

ii 

ii 

11-16 

H 

ii 

M 

ii 

11-18 

•i 

ii 

ii 

ii 

11-20 

ii 

ii 

ii 

ii 

12-12 

304. 

76 

203.78 

1.22 

55.37 

12-14 

ti 

ii 

ii 

ti 

12-16 

ii 

ii 

H 

ii 

12-18 

it 

H 

n 

ii 

12-20 

ii 

ii 

ii 

ii 

13-12 

269. 

11 

190.57 

1.15 

48.42 

13-14 

ii 

ii 

ii 

H 

13-16 

ii 

n 

ii 

ii 

13-18 

ii 

ii 

ii 

H 

13-20 

ii 

ii 

ii 

■i 

14-12 

234. 

36 

176.81 

1.08 

41.68 

14-14 

ii 

ii 

ii 

ii 

14-16 

ii 

H 

H 

■I 

14-18 

ii 

ii 

ii 

ii 

14-20 

n 

ii 

ii 

H 

-83- 


Table  A. 1.7  -  Dimensions  of  Member  #7 
(Factor  Safety  =  1.5) 


Frame 

(cm^ 

h 

t 

AF9 
(cmz ) 

(S2   -   S5) 

!) 

(cm) 

(cm) 

10-12 

372. 

38 

226. 

88 

1.33 

68.63 

10-14 

■I 

n 

n 

ii 

10-16 

n 

n 

ii 

ii 

10-18 

ii 

ii 

ii 

ii 

10-20 

ii 

ii 

ii 

ip 

11-12 

414. 

01 

240. 

07 

1.40 

76.84 

11-14 

ii 

ii 

ii 

" 

11-16 

ii 

n 

H 

ii 

11-18 

M 

ii 

n 

n 

11-20 

ii 

n 

H 

it 

12-12 

454. 

45 

252. 

25 

1.46 

84.84 

12-14 

n 

ii 

ii 

ii 

12-16 

ii 

M 

ii 

ii 

12-18 

ii 

ii 

H 

■I 

12-20 

H 

■I 

ii 

ii 

13-12 

494. 

78 

263. 

88 

1.52 

92.84 

13-14 

ii 

ii 

H 

n 

13-16 

ii 

ii 

ii 

n 

13-18 

ii 

it 

ii 

■1 

13-20 

H 

ii 

H 

ii 

14-12 

535. 

00 

275. 

01 

1.58 

100.84 

14-14 

ii 

ii 

H 

ii 

14-16 

ii 

ii 

ii 

ii 

14-18 

ii 

H 

ii 

ii 

14-20 

ii 

H 

ii 

ii 

-84- 


Table   A. 1.8    -   Dimensions   of  Member    #8 
(Factor   Safety  =    1.5) 


Frame 

T-. 

h 

t 

(cm2) 

(S2    -   S5) 

(cm2) 

(cm) 

(cm) 

10-12 

272.67 

188.33 

1.18 

49.26 

10-14 

221.67 

168.36 

1.08 

39.37 

10-16 

173.45 

147.23 

.97 

30.11 

10-18 

127.87 

124.34 

.85 

21.4  7 

10-20 

85.85 

99.26 

.72 

13.68 

11-12 

272.67 

188.33 

1.18 

49.26 

11-14 

221.67 

168.36 

1.08 

39.37 

11-16 

173.45 

147.23 

.97 

30.11 

11-18 

127.87 

124.34 

.85 

21.47 

11-20 

85.85 

99.26 

.72 

13.68 

12-12 

272.67 

188.33 

1.18 

49.26 

12-14 

221.67 

168.36 

1.08 

39.37 

12-16 

173.45 

147.23 

.97 

30.11 

12-18 

127.87 

124.34 

.85 

21.47 

12-20 

85.85 

99.26 

.72 

13.68 

13-12 

272.67 

188.33 

1.18 

49.26 

13-14 

221.67 

168.36 

1.08 

39.37 

13-16 

173.45 

147.23 

.97 

30.11 

13-18 

127.87 

124.34 

.85 

21.47 

13-20 

85.85 

99.26 

.72 

13.68 

14-12 

272.67 

188.33 

1.18 

49.26 

14-14 

221.67 

168.36 

1.08 

39.37 

14-16 

173.45 

147.23 

.97 

30.11 

14-18 

127.87 

124.34 

.85 

21.47 

14-20 

85.85 

99.26 

.72 

13.68 

-85- 


Table  A. 1.9  -  Dimensions  of  Member  #9 
(Factor  Safety  =  1.5) 


Frame 

AT-> 
(cm2) 

h 

t 

AFo 
(cmz) 

(S2  -  S5) 

(cm) 

(cm) 

10-12 

1278.87 

454.35 

2.47 

154.83 

10-14 

1031.25 

406.26 

2.23 

123.78 

10-16 

800.16 

355.82 

1.98 

94.96 

10-18 

581.14 

300.72 

1.70 

67.83 

10-20 

381.06 

240.29 

1.40 

43.30 

11-12 

1278.87 

454.35 

2.47 

154.83 

11-14 

1031.25 

406.26 

2.23 

123.78 

11-16 

800.16 

355.82 

1.98 

94.96 

11-18 

581.14 

300.72 

1.70 

67.83 

11-20 

381.06 

240.29 

1.40 

43.30 

12-12 

1278.87 

454.35 

2.47 

154.83 

12-14 

1031.25 

406.26 

2.23 

123.78 

12-16 

800.16 

355.82 

1.98 

94.96 

12-18 

581.14 

300.72 

1.70 

67.83 

12-20 

381.06 

240.29 

1.40 

43.30 

13-12 

1278.87 

454.35 

2.47 

154.83 

13-14 

1031.25 

406.26 

2.23 

123.78 

13-16 

800.16 

355.82 

1.98 

94.96 

13-18 

581.14 

300.72 

1.70 

67.83 

13-20 

381.06 

240.29 

1.40 

43.30 

14-12 

1278.87 

454.35 

2.47 

154.83 

14-14 

1031.25 

406.26 

2.23 

123.78 

14-16 

800.16 

355.82 

1.98 

94.96 

14-18 

581.14 

300.72 

1.70 

67.83 

14-20 

381.06 

240.29 

1.40 

43.30 

-86- 


Table  A. 1.10  -  Dimensions  of  Member  #1 
(Factor  Safety  =  2.0) 


Frame 

AT0 

h 

t 

(cm2) 

(S2   -   S5) 

(cm2) 

(cm) 

(cm) 

10-12 

514.90 

269.50 

1.55 

96.84 

10-14 

ii 

n 

ii 

ii 

10-16 

■i 

■I 

■1 

n 

10-18 

■I 

ii 

■i 

n 

10-20 

■I 

ii 

n 

ii 

11-12 

465.07 

255.36 

1.48 

86.95 

11-14 

■I 

ii 

ii 

ii 

11-16 

u 

ii 

H 

ii 

11-18 

n 

M 

H 

n 

11-20 

it 

■I 

n 

■1 

12-12 

417.21 

241.05 

1.41 

77.47 

12-14 

ii 

ii 

ii 

ii 

12-16 

u 

H 

■1 

n 

12-18 

ii 

ii 

ii 

ii 

12-20 

ii 

ii 

ii 

ii 

13-12 

372.38 

226.88 

1.33 

68.63 

13-14 

ii 

ii 

ii 

ii 

13-16 

■i 

ii 

ii 

n 

13-18 

ii 

ii 

H 

ii 

13-20 

ii 

ii 

ii 

ii 

14-12 

370.24 

226.18 

1.33 

68.21 

14-14 

364.89 

224.43 

1.32 

67.16 

14-16 

329.50 

212.50 

1.26 

60.21 

14-18 

n 

ii 

ii 

11 

14-20 

n 

H 

H 

H 

-87- 


Table  A. 1.11  -  Dimensions  of  Member  #2 
(Factor  Safety  =2.0) 


Frame 

AT 

h 

t 

AFo 

(S2  -  S5) 

(cm2) 

(cm) 

(cm) 

(cm2) 

10-12 

736.05 

325.25 

1.83 

141.05 

10-14 

770.63 

333.17 

1.87 

148.00 

10-16 

804.13 

340.66 

1.90 

154.74 

10-18 

840.74 

348.68 

1.94 

162.11 

10-20 

876.26 

356.30 

1.98 

169.26 

11-12 

790.53 

337.64 

1.89 

152.00 

11-14 

826.10 

345.50 

1.93 

159.16 

11-16 

862.68 

353.40 

1.97 

166.53 

11-18 

900.27 

361.36 

2.01 

174.11 

11-20 

939.91 

369.57 

2.05 

182.11 

12-12 

842.83 

349.13 

1.95 

162.53 

12-14 

881.48 

357.40 

1.99 

170.32 

12-16 

921.14 

365.70 

2.03 

178.32 

12-18 

960.76 

373.81 

2.07 

186.32 

12-20 

1001.39 

381.96 

2.11 

194.53 

13-12 

895.05 

360.26 

2.00 

173.05 

13-14 

936.78 

368.92 

2.04 

181.47 

13-16 

977.44 

377.18 

2.09 

189.68 

13-18 

1020.13 

385.66 

2.13 

198.32 

13-20 

1038.87 

389.33 

2.15 

202.11 

14-12 

927.40 

366.99 

2.03 

179.58 

14-14 

975.35 

376.76 

2.08 

189.26 

14-16 

1033.66 

388.32 

2.14 

201.05 

14-18 

1078.39 

396.96 

2.18 

210.11 

14-20 

1122.05 

405.23 

2.23 

218.95 

-88- 


Table  A. 1.12  -  Dimensions  of  Member  #3 
(Factor  Safety  =2.0) 


Frame 

Arp, 

h 

t 

AF 

(cm2) 

(S2   -   S5) 

(cm^ 

!) 

(cm) 

(cm) 

10-12 

511. 

11 

285.55 

1.53 

72.48 

10-14 

ii 

n 

ii 

ii 

10-16 

n 

■i 

ii 

ii 

10-18 

M 

n 

ii 

ii 

10-20 

■I 

ii 

ii 

n 

11-12 

462. 

13 

270.67 

1.46 

65.12 

11-14 

ii 

ii 

ii 

ii 

11-16 

H 

H 

ii 

ii 

11-18 

ii 

■i 

it 

ii 

11-20 

ii 

H 

ii 

H 

12-12 

415. 

12 

255.62 

1.39 

58.08 

12-14 

ii 

ii 

ii 

it 

12-16 

ii 

■■ 

■I 

n 

12-18 

ii 

ii 

n 

ii 

12-20 

n 

ii 

H 

ii 

13-12 

370. 

07 

240.37 

1.32 

51.36 

13-14 

ti 

ii 

ii 

ii 

13-16 

H 

it 

ii 

ii 

13-18 

it 

H 

ii 

ii 

13-20 

ii 

n 

H 

u 

14-12 

134. 

63 

138.13 

.84 

16.96 

14-14 

161. 

59 

152.97 

.91 

20.80 

14-16 

326. 

98 

224.90 

1.25 

44.96 

14-18 

ii 

ii 

ii 

ii 

14-20 

n 

ii 

n 

ii 

-89- 


Table  A. 1.13  -  Dimensions  of  Member  #5 
(Factor  Safety  =  2.0) 


Frame 

(cm2) 

h 

t 

AF 
(cm2) 

(S2  -  S5) 

(cm) 

(cm) 

10-12 

1149.83 

449.13 

2.22 

151.29 

10-14 

1462.79 

508.97 

2.49 

194.29 

10-16 

1792.07 

565.35 

2.74 

239.71 

10-18 

2126.50 

617.52 

2.98 

286.00 

10-20 

2467.32 

666.62 

3.20 

333.29 

11-12 

1133.80 

444.00 

2.11 

134.86 

11-14 

1488.72 

513.62 

2.51 

197.86 

11-16 

1815.84 

569.21 

2.76 

243.00 

11-18 

2148.15 

620.75 

2.99 

289.00 

11-20 

2486.87 

669.33 

3.21 

336.00 

12-12 

1206.09 

460.43 

2.27 

159.00 

12-14 

1512.57 

517.87 

2.53 

201.14 

12-16 

1837.54 

572.71 

2.78 

246.00 

12-18 

2167.73 

623.66 

3.01 

291.71 

12-20 

2505.38 

671.88 

3.22 

338.57 

13-12 

1230.03 

465.17 

2.29 

162.29 

13-14 

1533.31 

521.54 

2.55 

204.00 

13-16 

1856.14 

575.70 

2.79 

248.57 

13-18 

2185.25 

626.25 

3.02 

294.14 

13-20 

2520.80 

674.01 

3.23 

340.71 

14-12 

1322.61 

483.05 

2.37 

175.00 

14-14 

1609.98 

534.88 

2.61 

214.57 

14-16 

1872.67 

578.34 

2.80 

250.86 

14-18 

2200.71 

628.53 

3.03 

296.29 

14-20 

2535.20 

675.98 

3.24 

342.71 

-90- 


Table  A. 1.14  -  Dimensions  of  Member  #6 
(Factor  Safety  =2.0) 


Frame 

"TV 

(cm' 

h 

t 

AF 
(cm^) 

(S2    -   S5) 

!) 

(cm) 

(cm) 

10-12 

459. 

76 

253. 

81 

1.47 

85.89 

10-14 

ii 

ii 

ii 

ii 

10-16 

ii 

H 

ii 

ii 

10-18 

H 

ii 

n 

n 

10-20 

n 

n 

ii 

n 

11-12 

410. 

82 

239. 

08 

1.40 

76.21 

11-14 

ii 

n 

ii 

ii 

11-16 

H 

ii 

" 

ii 

11-18 

ii 

n 

ii 

n 

11-20 

ii 

H 

■1 

ii 

12-12 

364. 

89 

224. 

43 

1.32 

67.16 

12-14 

ii 

ii 

ii 

ii 

12-16 

■■ 

ii 

n 

ii 

12-18 

ii 

ii 

" 

n 

12-20 

ii 

ii 

ii 

ii 

13-12 

320. 

90 

209. 

51 

1.25 

58.53 

13-14 

n 

ii 

ii 

ii 

13-16 

ii 

ii 

ii 

ii 

13-18 

ii 

H 

ii 

H 

13-20 

H 

H 

ii 

n 

14-12 

279. 

93 

194. 

67 

1.17 

50.53 

14-14 

ii 

ii 

ii 

ii 

14-16 

ii 

ii 

ii 

ii 

14-18 

ii 

n 

ii 

n 

14-20 

ii 

ii 

ii 

H 

-91- 


Table  A. 1.15  -  Dimensions  of  Member  #7 
(Factor  Safety  =2.0) 


Frame 

(cm2) 

h 

t 

AF 
(cm2) 

(S2   -   S5) 

(cm) 

(cm) 

10-12 

445.95 

249.74 

1.45 

83.16 

10-14 

ii 

n 

11 

H 

10-16 

ii 

ii 

ti 

M 

10-18 

H 

ii 

II 

ii 

10-20 

ii 

ii 

II 

ir 

11-12 

495.84 

264.18 

1.52 

93.05 

11-14 

ii 

n 

ii 

M 

11-16 

n 

ii 

■I 

ii 

11-18 

ii 

ii 

ii 

n 

11-20 

ii 

ii 

n 

M 

12-12 

545.57 

277.87 

1.59 

102.95 

12-14 

■i 

ii 

n 

ii 

12-16 

n 

ii 

ii 

■I 

12-18 

M 

n 

ii 

n 

12-20 

ii 

ii 

ii 

H 

13-12 

594.13 

290.64 

1.65 

112.63 

13-14 

it 

ii 

ii 

ii 

13-16 

ii 

■I 

ii 

ii 

13-18 

ii 

ii 

H 

n 

13-20 

ii 

ii 

H 

ii 

14-12 

642.58 

302.88 

1.71 

122.32 

14-14 

ti 

ii 

ii 

ii 

14-16 

ii 

H 

ii 

ii 

14-18 

ii 

ii 

n 

ii 

14-20 

n 

ii 

ii 

ii 
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Table  A. 1.16  -  Dimensions  of  Member  #8 
(Factor  Safety  =2.0) 


Frame 

(cm2) 

h 

t 

AF, 

(S2  -  S5) 

(cm) 

(cm) 

(cm2) 

10-12 

326.50 

207.48 

1.28 

59.79 

10-14 

265.10 

185.50 

1.17 

47.79 

10-16 

206.39 

161.94 

1.04 

36.42 

10-18 

152.42 

137.10 

.91 

26.11 

10-20 

101.89 

109.43 

.77 

16.63 

11-12 

326.50 

207.48 

1.28 

59.79 

11-14 

265.10 

185.50 

1.17 

47.79 

11-16 

206.39 

161.94 

1.04 

36.42 

11-18 

152.42 

137.10 

.91 

26.11 

11-20 

101.89 

109.43 

.77 

16.63 

12-12 

326.50 

207.48 

1.28 

59.79 

12-14 

265.10 

185.50 

1.17 

47.79 

12-16 

206.39 

161.94 

1.04 

36.42 

12-18 

152.42 

137.10 

.91 

26.11 

12-20 

101.89 

109.43 

.77 

16.63 

13-12 

326.50 

207.48 

1.28 

59.79 

13-14 

265.10 

185.50 

1.17 

47.79 

13-16 

206.39 

161.94 

1.04 

36.42 

13-18 

152.42 

137.10 

.91 

26.11 

13-20 

101.89 

109.43 

.77 

16.63 

14-12 

326.50 

207.48 

1.28 

59.79 

14-14 

265.10 

185.50 

1.17 

47.79 

14-16 

206.39 

161.94 

1.04 

36.42 

14-18 

152.42 

137.10 

.91 

26.11 

14-20 

101.89 

109.43 

.77 

16.63 
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Table  A. 1.17  -  Dimensions  of  Member  #9 
(Factor  Safety  =2.0) 


Frame 

™r„ 

h 

t 

AF 
(cm2) 

(S2  -  S5) 

(cm2) 

(cm) 

(cm) 

10-12 

1539.02 

500.09 

2.70 

187.57 

10-14 

1240.44 

447.21 

2.44 

150.00 

10-16 

960.29 

391.43 

2.16 

114.91 

10-18 

698.25 

331.27 

1.86 

82.30 

10-20 

455.88 

264.41 

1.52 

52.43 

11-12 

1539.02 

500.09 

2.70 

187.57 

11-14 

1240.44 

447.21 

2.44 

150.00 

11-16 

960.29 

391.43 

2.16 

114.91 

11-18 

698.25 

331.27 

1.86 

82.30 

11-20 

455.88 

264.41 

1.52 

52.43 

12-12 

1539.02 

500.09 

2.70 

187.57 

12-14 

1240.44 

447.21 

2.44 

150.00 

12-16 

960.29 

391.43 

2.16 

114.91 

12-18 

698.25 

331.27 

1.86 

82.30 

12-20 

455.88 

264.41 

1.52 

52.43 

13-12 

1539.02 

500.09 

2.70 

187.57 

13-14 

1240.44 

447.21 

2.44 

150.00 

13-16 

960.29 

391.43 

2.16 

114.91 

13-18 

698.25 

331.27 

1.86 

82.30 

13-20 

455.88 

264.41 

1.52 

52.43 

14-12 

1539.02 

500.09 

2.70 

187.57 

14-14 

1240.44 

447.21 

2.44 

150.00 

14-16 

960.29 

391.43 

2.16 

114.91 

14-18 

698.25 

331.27 

1.86 

82.30 

14-20 

455.88 

264.41 

1.52 

52.43 
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Table  A. 1.18 
Limit  Design  Frame  Weights  for  Single  Cross-Tie  Tanker 


Frame  Weight         Frame  Weight 
Frame  (metric  tons)         (metric  tons) 

(S2  -  S5)  FS  =  1.5  FS  =  2.0 


10-12 
10-14 
10-16 
10-18 
10-20 
11-12 
11-14 
11-16 
11-18 
11-2  0 
12-12 
12-14 
12-16 
12-18 
12-20 
13-12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


41. 

,77 

42. 

,71 

44. 

,82 

48. 

,96 

56. 

,08 

41. 

,55 

41. 

,99 

44. 

,41 

48. 

,67 

54. 

,77 

41. 

,47 

41. 

,90 

44. 

,29 

48. 

,53 

54. 

.57 

41. 

,72 

42. 

,14 

44, 

,50 

48. 

,68 

54, 

.73 

41, 

.67 

42, 

.89 

44. 

,92 

49, 

.11 

55, 

.09 

48, 

,89 

49, 

,09 

51, 

,77 

56, 

.44 

63, 

,45 

47, 

.85 

48, 

.63 

51, 

.13 

55, 

,94 

62, 

,77 

48, 

.41 

48, 

,50 

50. 

,96 

55. 

,59 

62, 

,48 

48, 

.65 

48, 

.70 

51, 

.09 

55, 

.71 

62, 

.34 

48, 

.51 

48, 

.40 

51, 

.54 

56, 

.13 

62, 

.85 
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APPENDIX  B  -  TANKER  FRAME  WITH  RIGID  CORNERS 

In  large  ships  which  have  relatively  simple 
web  frame  geometries,  the  webs  of  the  individual  beams 
are  likely  to  be  quite  deep.   Brackets  are  also  rather 
large. 

In  this  study,  one  such  frame  is  examined 
using  the  assumption  of  rigid  extensions  in  the  bracket 
and  beam  overlap  regions.   The  web  frame  numbered 
OC22,49.2t  was  selected  from  Lund's  paper  [5]  as  a 
representative  frame  from  which  a  model  could  be 
developed.   This  frame  was  among  the  lighter  frames 
having  zero  cross-ties  and  is  shown  in  Figure  10  of 
Reference  f5]. 

The  frame  model  is  shown  in  Figure  A. 3.   The 
dots  on  the  system  lines  represent  the  end  points  of 
the  rigid  extensions.   Here,  as  in  Appendix  A.,  the 
system  lines  are  assumed  to  be  defined  on  the  inside 
surface  of  the  shell  plating  and  on  the  center  tank 
face  of  the  longitudinal  bulkhead  plating.   Plastic 
neutral  axes  of  the  individual  beam  elements  calculated 
did  in  fact  lie  on  or  very  near  these  lines.   The 
principal  dimensions  of  the  ship  under  consideration  are 
the  same  as  those  listed  in  Appendix  A. 

The  loading  conditions  considered  are  those 
defined  by  Lund  [5],  discussed  in  Appendix  A  and  shown 
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in  Figure  A. 2. 

The  assumed  lengths  of  bracket  legs  are 
identical  to  those  used  by  Lund  and  are  listed  in  Table 
2  of  Reference  [5]. 

The  lengths  of  the  rigid  extensions  were 
calculated  by  hand  using  the  procedures  outlined  in 
Section  2.2  and  Reference  [6].   Although  Lund  only 
calculates  these  values  once  for  each  frame,  he  does 
not  list  his  results.   It  is  therefore  assumed,  but 
not  known  as  fact  that  rigid  extension  lengths  used  in  this 
study  are  approximately  equal  to  values  used  by  Lund  for 
the  same  frame.   The  actual  values  used  are  listed  in 
Table  A. 2.1. 

In  this  analysis  the  calculations  required 
to  develop  the  limit  design  problem  were  done  by  hand. 
The  assumption  of  rigid  joint  regions  required  that 
each  proposed  collapse  mechanism  be  carefully  inspected 
in  insure  kinematic  admissibility.   Joint  rotations 
were  not  considered  as  fundamental  mechanisms  in  the 
formations  of  mechanisms  which  involved  more  than  one 
frame  element.   A  discussion  of  the  joint  rotation 
problem  is  given  in  Section  2.3. 

The  results  of  the  hand  calculations  are 
listed  in  Table  A. 2. 2   The  values  of  the  constraint 
matrix  coefficients  (i.e.  internal  energy  dissipation 
coefficients)  as  well  as  the  values  of  the  constraint 
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equation  right  hand  sides  (i.e.  external  work  rates) 
are  tabulated. 

It  should  be  noted  that  although  the  values 
of  constraint  relations  are  shown  in  the  table  to  be 
segregated  by  load  condition,  they  were  combined  as  one 
system  in  the  linear  programming  problem.   This  allowed 
the  system  to  be  truly  optimized  for  all  considered 
loadings.   In  general,  this  is  not  possible  in  the 
elastic  design  methods.   Most  elastic  procedures  optimize 
for  each  load  condition  and  then  select  the  larger 
dimensions  required  for  each  element. 

The  linear  programming  subroutine  M.P.S.X. 
is  used  to  optimize  the  results  of  the  hand  calculations 
as  discussed  in  Appendix  C.   Again,  as  in  Appendix  A, 
the  output  from  the  procedure  is  a  listing  of  minimum 
values  of  plastic  limit  moment  required  to  make  the 
structure  safe  under  the  prescribed  loading  conditions. 

The  dimensions  of  the  beam  elements  were 
calculated  based  on  the  relation 


where 


Mo  =  ay  zp 


z   =  A  a 


In  these  relationships  M0  is  the  plastic  limit 
moment  calculated  in  the  limit  design  problem,  a      is  the 
yield  stress  of  the  material,  A  is  the  total  required  area 
of  the  cross  section  and  a  is  the  distance  between 
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centroids  of  the  two  equal  areas  lying  above  and  below  the 
plastic  neutral  axis. 

In  order  to  allow  comparison  to  the  parent 
Lund  frame,  the  beam  elements  are  sized  to  meet  the 
plastic  moment  requirements  based  upon  the  proportions 
of  the  beams  which  Lund  used  as  starting  values  for 
his  frame  iterations.   For  this  frame,  starting  values  of 
web  depth,  web  thickness  and  top  flange  area  are  shown 
in  Table  A. 2. 3.   Also  shown  are  values  of  plate  area 
which  are  assumed  to  remain  constant. 

Values  of  beam  geometry  were  calculated  using 
factors  of  safety  of  1.0,  1.5,  2.0  and  2.5.   A  corrosion 
allowance  of  0.2  cm  was  applied.   Table  A. 2. 4  gives 
resulting  values  of  net  beam  cross  sectional  area  A  , 
web  depth  h,  web  thickness  t  and  top  flange  area  AF. 

Frame  weights  were  calculated  for  one  half  of 

complete  frames  for  each  of  the  factors  of  safety 

considered.   The  relation 

6 
W  -  p  I    {1±  hj,   -  ^5  (C.  +  ck)AT} 

was  used  to  calculate  these  weight  values  where  p  is 

material  density,  1.  is  the  member  length  and  c.  and  c, 

*'      i  3  k 

are  the  overlaps  at  the  respective  beam  ends.   This 
expression  was  also  used  by  Lund.   Consistent  with  Lund's 
calculations,  weight  of  web  stiff eners  and  brackets  are 
neglected. 
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Values  of  frame  weight  are  plotted  against 
factors  of  safety  for  this  analysis  and  are  shown  in 
Figure  A. 6.   Also  plotted  is  the  weight  of  the  Lund 
elastic  design  frame  at  a  factor  of  safety  of  1.5. 
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Table  A. 2.1  -  Rigid  Extension  Lengths 


Rigid 
Extension 

Length 
(cm) 

LB 

504 

Lc 

630 

LD 

128 

LE 

428 

L„ 

553 

Rigid 
Extension 

Length 
(cm) 

LM 

453 

Ln 

453 

Lo 

153 

LP 

392 

L^ 

382 

-101- 


x 
w 


oor~-r^r-oooor-ooooc^-r-~ooooooix> 


Oi 

c 
o 

•H 

-p 

D1 
ft) 

c 
On 

•H 

to 

P 

+J 

■H 

e 

•H 

IH 

O 

tO 
4J 
C 

Q) 

•H 

u 

■H 

m 
O 

u 


CO 


s 


i-H'3'Or~cNcoinrr><7\r-HcriCN(N.— ir-- 
nr^rviincMrHr-iHiHrHinrominnH 


in 

in 


IT) 

s 


o 

<D  00 

TJ<  "<3<  CM 

ic 

r-  r- 

Nroin 

o 

r^  <* 

(N0OH 

<j\ 

O  00 

WOrt 

■T  CN 


CN 


CN 


s 


in 

*r  ^r 

oo 

r»  m 

en 

cti  <n 

r- 

r-  r» 

CM 

o 

m 

CN 

<D 

ai 

vo 

in 

m 


ro 


vo 

r- 

CM 

O   <£> 

00 

** 

<* 

H 

V£>   CTi 

r^- 

m 

00 

rH 

o  cr\ 

^r 

r- 

r- 

-H 

<J>  00 

00 

♦ 

• 

• 

•       • 
CN  CN 

• 

CN 
< 

•9 

Eh 


(0 

6 
to 

•H 
C 

rd 

& 


CN 

s 


O       00 

co  o  o  r^  o 


W5 

en 
oo 


m 

o 
ro 


CN 


cd 
o 


o  o 


00  <X> 

cn  in 

CN  -^ 
^  00 


m 
oo 


rfJcQCjaWfaUffiHinK^SSO 


-102- 


w 


^cococncocooocooor^oooooooooo 


O'*<N(Ti<Nr~Tj'T}''rrOVD(TiLn00U3 

i,(Nwooi^ovooooor«0'*oair^ 


VD 


ro 

in 

o  o        LD 


"*  iH 


lo 


r- 

<J< 

«# 

CO 

o 

ro 

CN 

if) 

ro 

00 

cn 

H 

U3 

ro 

ro 

H 

CM 


S 


•^r 

^r 

r~ 

r- 

Ch 

Ol 

r> 

p» 

CN  CN  O  M" 

rH  ro  CN  <-\ 

00  >>0  CTi  CN 

•    •  •    • 

rH  ro  ro  CN 


a) 

3 
c 

■H 

-P 

C 

o 
u 


ro 


n  r-  id  vo  r~~  oo 

Cfl  WKl'JOO 

HCTi^rnm  in 

<J\  CM  00  M>  l£>  <N 
rH  n  CN  m 


CN 
CN 


rH 

(0 

En 


CO 

a 

w 

■H 

c 

A 
o 

i 

CN 

rrj 


CN 


^       ro  r» 

o      rH  en 

nocniNO 

rH  ^  cn  cn  ro 


CO 

o 
in 

CN 


o 
ro  o 

•     • 
•-{   CN 


m 


o  o 


CO 


"^  rH 


<«UQWhUKHh«rlS20 


-103- 


T3 
<D 

a 

•H 
4J 

C 

o 
u 


(N 
ON 

CD 

rH 

Eh 


A, 
X 

w 


1X1 


S 


in 


ro 

g 


cor-r^oooo^rroooocnoooooooooooo 


inoonvDvoin[^ix>cHroo>ro<Tir~-in(N 
vocNOOooor-Hoco^rcorom'tfooo 


00 

ro  o 

r^ 

O  rH 

ro 

<N  (N 

rH 

O   O 

VD   O 

o 

CO  VD 

H 

rH  M 

CN 

O  rH 

o 

O  in  CN 

oo  oo  ro 

^  <Ti  r^ 

•    •    • 

rr  vd  cn 


00  ^D  <J3  «* 

cn  o  o  ro 

nji  ^}i  «g<  CA 

o  o  o  r- 

•      •       •  • 

■H  "3<  ^  iH 


CO^fOO 

oo  ro 

00 

in  oo  ^  i-» 

CN  o 

in 

in  oo  vo  m 

■sr  cn 

^r> 

tOOOH 

o  o 

CN 

CN  CN  ^r 


CN  CN 


w 

e 

en 

■** 

00 

-H 

ro 

in 

§ 

CN 

o 

m 

g 

ro 

o  o  00 

X 

• 

•    •    • 

o 

<-i 

■»j<  ro  cn 

0 

g 

ro 

=tfc 

VI5 

in 

13 

rH 

»* 

CO 

g 

o  o 

00 

O 

•       • 

• 

3 

**  f-\ 

f-f 

<-\ 

00 

r- 

ro 

ro 

c^ 

<-\ 

in 

o  o 

O  O 

■ 

« 

•       • 

•      • 

CN 

00 

CN    r-t 

rH    CN 

o 

CD 


CN 

rH 


o  o 

rH    r-\ 


<fflUQWhOXHo«J220P< 


-104- 


CM 

w 


oo  oo  oo  oo  vo 


en 


H  ID(NCC(N 

oo  in  oo  vo  vo 

O  O  ID  CTi  <T\ 

HHinnin 


VD 


m 


m 


T3 
0) 

C 

•H 
-P 

c 
o 
u 


r-\ 

(0 
Eh 


g 

CO 
•H 

c 

u 

a 
■a 

n5 
0 
3 


s 


<N 


o 
r- 


ro 
O 
ro  O 

•     • 


o  o 

•       • 

■^  H 


<  m  u  Q  w 


-105- 


Table  A. 2. 3  -  Initial  Member  Dimensions 


Web 

Top  Flange 

Plate  Flange 

Member 

Height 

Thickness 

Area 

Area 

No. 

(cm) 

(cm) 

(cm2) 

(cm2) 

1 

400 

2.0 

120 

750 

2 

400 

1.8 

120 

750 

3 

300 

1.5 

120 

660 

4 

300 

1.5 

120 

420 

5 

300 

1.4 

80 

750 

6 

240 

1.25 

80 

750 
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Table  A. 2. 4  -  Dimensions  of  Limit  Design  Results 


AT 

321. 

00 

1130. 

.00 

415. 

.00 

533, 

,00 

40. 

00 

98. 

00 

475. 

,02 

1445. 

,16 

601. 

,51 

640. 

,47 

71. 

,16 

123. 

,38 

596. 

,51 

1739. 

,34 

722. 

,95 

769. 

,03 

84. 

,32 

153. 

,95 

FS  =  1.0 

1  321.00  236.52  1.28  42.07 

2  1130.00  463.94  2.29  161.43 

3  415.00  255.64  1.48  87.35 

4  533.00  264.56  1.52  93.03 

5  40.00  84.77  .60  6.39 

6  98.00  96.01  1.09  12.65 

FS  =  1.5 

1  475.02  270.27  1.55  54.78 

2  1445.16  505.78  2.48  191.86 

3  601.51  292.54  1.66  114.11 

4  640.47  302.36  1.71  121.89 

5  71.16  95.81  .65  8.16 

6  123.38  121.88  .83  20.63 

FS  =  2.0 

1  596.51  296.53  1.69  66.39 

2  1739.34  556.69  2.71  232.43 

3  722.95  322.21  1.81  138.42 

4  769.03  332.80  1.86  147.68 

5  84.32  105.78  .69  9.95 

6  153.95  137.86  .92  26.39 

FS  =  2.5 

1  655.07  320.00  1.80  76.96 

2  2009.29  599.74  2.90  269.76 

3  833.37  347.08  1.94  160.62 

4  887.19  358.61  1.99  171.47 

5  96.14  114.09  .73  11.57 

6  176.69  148.74  .97  30.73 


1. 

,28 

2. 

.29 

1. 

.48 

1. 

.52 

.60 

1, 

.09 

1, 

.55 

2, 

,48 

1, 

,66 

1 

.71 

.65 

i 

.83 

1, 

.69 

2, 

.71 

1 

.81 

1, 

.86 

i 

.69 

< 

.92 

1, 

.80 

2, 

.90 

1, 

.94 

1, 

.99 

.73 

i 

.97 

AF 
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APPENDIX  C  -  NUMERICAL  ROUTINE 

C.l   Use  of  the  M.P.S.X. 

The  M.P.S.X.  linear  programming  subroutine  used 
in  this  analysis  was  found  to  be  very  sensitive  to  format 
errors.   The  user's  guides  and  instructions  available  on 
the  system  are  not  very  careful  to  point  out  possible 
pit  falls.   It  is  the  purpose  of  this  section  to  list 
the  control  language  statements  used  in  solving  the 
linear  programming  problems  of  this  analysis  and  to 
point  out  general  problems  that  were  encountered  in 
using  the  M.P.S.X..   Before  attempting  to  use  this 
subroutine  the  user  should  become  very  familiar  with 
References  [31  -  33]. 

A  sample  listing  of  the  control  program  with 

sample  data  points  follows: 

//LP   EXEC  MPSX 
//MPSCOMP.SYSIN   DD   * 

PROGRAM 

INITIALZ 

MOVE ( XDATA , ' LUND ' ) 

MOVE (XPBNAME , ' PBFILE ' ) 

CONVERT ( ' SUMMARY ' ) 

BCDOUT 

SETUP 

MOVE(XOBJ, 'MINWT* ) 

MOVE ( XRHS , ' WORK ■ ) 

CRASH 

PRIMAL 

SOLUTION 

RANGE 

EXIT 

PEND 

/* 

//MPSEXEC.SYSIN   DD   * 
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NAME 
ROWS 
G   MECHA 

COLUMNS 
MOM1 


LUND 


RHS 


WORK 


MECHA 


MECHA 


4.0 
3.605E+8 


MECHB 


MECHB 


1.0 


4.228E+7 


ENDATA 
/* 


Explanations  of  the  control  language  statements 
can  be  found  in  Reference  [33].   The  format  sensitivity 
of  the  M.P.S.X.  has  been  mentioned  previously.   The 
greatest  difficulty  occurs  when  entering  data.   Page 
2  09  of  Reference  [33]  will  be  found  to  be  invaluable 
to  any  prospective  system  user.   This  page  tabulates  all 
required  data  entries  in  the  positions  required  by  the 
M.P.S.X..   The  data  set  name  and  all  data  vector  names 
must  be  left  justified  in  their  fields.   Failure  to 
have  the  data  set  name  left  justified  will  result  in 
having  the  data  set  not  read.   If  a  row  or  vector  name 
is  not  left  justified,  its  associated  data  point  will 
not  be  read.   If  a  row  or  vector  name  extends  beyond  the 
left  bound  of  its  field,  the  data  will  be  read  and  labeled 
with  an  erroneous  name.   Any  of  these  errors  associated 
with  row  and  vector  names  will  result  in  erroneous 
solutions. 

The  scale  values  of  the  data  may  be  placed 
anywhere  in  the  prescribed  field. 
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Finally,  Chapter  4  of  Reference  [33]  should  be 
checked  to  insure  that  data  and  vector  names  use  only 
"legal"  characters.   The  use  of  blanks  and  certain  other 
characters  will  result  in  abandoned  runs. 
C.2   M.P.S.X.  Costs 

The  two  positive  advantages  of  the  M.P.S.X. 
are  its  power  and  efficiency.   The  linear  programming 
problems  solved  in  this  analysis  are  moderately  large, 
having  roughly  fifty  constraint  relations  and  nine 
variables.   There  was  originally  some  fear  that  the 
solutions  of  these  problems  would  be  quite  expensive . 

It  was  found  that  problems  having  approximately 
5  0  constraints  and  6  variables  required  approximately 
0.325  minutes  in  total  job  run  time  at  a  cost  of 
approximately  four  dollars.   These  numbers  are  based 
on  computation  performed  in  solving  the  linear  programming 
problem  of  Appendix  B. 

Run  time  and  associated  costs  were  found  to  be 
very  closely  related  to  the  total  number  of  constraints. 
Changes  in  the  number  of  solution  variables  had  little  or 
no  effect  on  computation  time. 

The  total  cost  of  running  the  limit  design 
program  in  Appendix  A  for  twenty-five  different  frames 
including  the  card  punching  cost  was  $17.57.   The  cost 
of  solving  the  programming  problems  formed  by  the  FORTRAN 
program  was  $74.08.   This  resulted  in  a  total  computation 
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cost  of  $91.65  for  limit  design  problem  solutions  of 
twenty-five  frame  configurations. 
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17U70 

TK455  S     K1«9hornt  of  ship 

Averse  web  frames. 


Thesis  171170 

K4d5  Kinghorn 

Limit   design   of   ship 

♦:ransverse  web   *r3rres. 


■design  of  ship  transverse  web  ham 


Limit  < 
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